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Calculators  for  R-2920-ARPA, 

Electromagnetic  Properties  and  Communication  Characteristics 

of  PACSAT 


Three  calculators  have  been  developed  to  supplement  the  plots 
and  calculations  presented  in  this  report. 

(1)  Passive  Satellite  Calculator.  This  calculator  is  a  circular 
slide  rule  that  permits  numerical  evaluations  of  typical  passive-satellite 
link  parameters.  It  solves  Eq.  (A. 25)  to  produce  results  such  as  those 
listed  in  Table  4.4. 

(2)  Satellite  Range  Calculator.  This  (and  the  following)  calculator 
uses  the  polar  stereographic  projection  of  the  Northern  Hemisphere  shown 
in  Fig.  4.5  to  calculate  the  slant  range  from  a  satellite  in  a  circular 
equatorial  orbit  to  an  arbitrary  point  within  its  field  of  view.  A  set 

of  overlays  covers  the  range  of  orbital  altitudes  from  3000  to  10,000  km. 
Unlike  Fig.  4.12,  which  is  based  on  relative  satellite  and  user  location, 
this  slide  rule  uses  actual  geographical  locations  and  facilities  visu¬ 
alizing  the  coverage  area. 

(3)  Satellite  Coverage  Calculator.  Using  the  same  projection  and 
satellite  orbits  as  the  Range  Calculator,  this  calculator  uses  a  pair  of 
overlays,  similar  to  Fig.  4.4,  to  depict  the  extremes  of  satellite  longi¬ 
tudes  for  mutual  visibility  between  a  specific  user  pair  or  for  users 
within  a  given  user  area.  Auxiliary  scales  yield  the  arc  of  mutual  visi¬ 
bility  and  the  mutual  viewing  time  for  a  generalized  pair  of  users 
(Figs.  4.7  and  4.8  give  these  quantities  for  specific  user  pairs). 

The  Passive  Satellite  Calculator  is  general  in  that  it  can  be  used 
over  a  wide  range  of  parameters  for  a  passive  satellite  system  with  no 
restrictions  on  the  orbit.  The  Satellite  Range  and  Coverage  Calculators 
are  specific  for  satellites,  whether  active  or  passive,  in  circular 
equatorial  orbits.  These  calculators  are  available  at  cost  from  The  Rand 
Corporation.  Inquiries  should  be  directed  to  the  Engineering  and  Applied 
Sciences  Department. 
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PREFACE 


This  report  presents  an  analysis  of  various  technical  and  system 
properties  of  a  passive  communication  satellite  concept  known  as 
PACSAT.  Included  are  electromagnetic  scattering,  orbital  and  geomet¬ 
rical  considerations,  and  communication  characteristics.  A  general 
description  of  PACSAT  and  its  potential  application  to  the  command  and 
control  of  strategic  forces  is  given  in  a  companion  Rand  Note:  Edward 
Bedrosian,  PACSAT:  A  Passive  Communication  Satellite  for  Survivable 
Command,  and  Control ,  N-1780-ARPA,  November  1981. 

The  research  reported  here  was  sponsored  by  the  Strategic  Tech¬ 
nology  Office,  Defense  Advanced  Research  Projects  Agency  (contract  MDA 
903-78-C-0281),  under  the  supervision  of  Dr.  Sherman  Karp.  It  is  an 
element  of  a  broad  range  of  strategic  communication  studies  being 
conducted  for  DARPA  by  Rand. 


SUMMARY 


This  report  presents  a  detailed  electromagnetic  and  communication 
system  analysis  of  a  passive  communication  satellite  known  as  PACSAT. 
The  satellite  is  a  linear  array  of  scatterers  oriented  toward  the 
center  of  the  earth  by  the  gravity-gradient  force.  The  scatterers  are 
assumed  to  consist  of  uniformly  spaced,  perfectly  conducting  spheres. 

First,  the  scattering  cross  section  of  the  array  is  developed  and 
the  effects  of  flexure  and  near-field  operation  are  considered.  The 
array  is  then  assumed  to  be  in  a  circular  equatorial  orbit  and  various 
system  and  orbital  parameters  are  calculated  for  users  in  a  service 
area  approximating  the  United  States  in  shape.  Finally,  communication 
characteristics  are  investigated  and  power  budgets  presented  for  typi¬ 
cal  system  configurations. 
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GLOSSARY 


The  defining  equation  (in  italics)  or  the  equation  of  first  ref¬ 
erence  is  given  at  the  right.  Figure  numbers  are  given  when  appro¬ 
priate.  Unit  vectors  are  formed  by  placing  bars  over  the  appropriate 
coordinate  designation. 

English 


A 

azimuth  of  line  of  sight 

Eq.  (4.6) 

Fig.  4.1 

Aeff 

effective  area 

Eq.  (2.35) 

An 

coefficient  of  nth  term  in  Mie 
series  expansion 

Eq.  (2.6) 

a 

sphere  radius 

Eq.  (2.5) 

axy 

scattering  coefficient,  y  to  x 
polarization 

Eq.  (2.11) 

Bn 

coefficient  of  nth  term  in  Mie 
expansion 

Eq.  (2.6) 

C(z) 

cosine  Fresnel  integral 

Eq.  (3.12) 

cn 

coefficient  of  nth  term  in  array 
length  expansion 

Eq.  (3.2) 

D 

maximum  deflection  of  flexed  array 

Eq.  (3.29) 

Fig.  3.4 

D 

diameter  of  parabolic  dish 

Eq.  (4.24) 

D 

directivity  pattern 

Eq.  (2.42) 

3 

normalized  flexure  parameter 

Eq.  (3.33) 

dk 

distance  from  array  center  to  kth 
element 

Eq.  (3.20) 

Fig.  3.3 

1*1 

magnitude  of  electric  field  vector 

Eq.  (2.19) 

E 

elevation  angle  of  line  of  sight 

Eq.  (4.8) 

Fig.  4.1 

Eb 

received  energy  per  bit 

Eq.  (4.19) 

incident  electric  field  vector 

Eq.  (2.1) 

E8 

scattered  electric  field  vector 

Eq.  (2.3) 

fl 

frequency  to  direct  first-order 
grating  lobe  toward  receiver 

Eq.  (2.57) 

f3dB 

3  dB  beamwidth  in  frequency 

Eq.  (2.68) 

G 

array  gain 

Eq.  (2.27) 
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antenna  gain 

orbital  altitude  (above  surface 
of  earth) 

geometrical  parameter 

spherical  Bessel  function  of  the 
third  kind 

defined  integral 

geometrical  parameter 

real  part  of  spherical  Bessel 
function  of  third  kind 

order  of  grating  lobe 

index  to  kth  array  element 

Boltzmann's  constant 
(1.38054  x  10~3J/°K) 

wavelength  parameter 

array  length 

link  loss  factor 

path  length  from  transmitter  to 
receiver  via  kth  element 

projected  array  length 

path  length  to  kth  element 

system  margin 

number  of  elements  In  array 

system  noise  spectral  density 

received  and  transmitted  powers 

associated  Legendre  function 

line-of-sight  path  lengths  to 
receiver  and  transmitter 

radius  of  the  earth  (6378  km) 

distance  to  reference  field  point 

radius  of  orbit 

sine  Fresnel  integral 

functions  in  Mie  series  expansion 

sine  integral 

array  element  spacing 

distance  along  array  to  kth 
element 


Eq.  (4.22) 
Eq.  (3.17) 

Eq.  (3.25) 
Eq.  (2.6) 

Eq.  (2.30) 
Eq.  (3.25) 
Eq.  (2.6) 

Eq.  (2.22) 
Eq.  (3.2) 
Eq.  (4.18) 

Eq.  (2.1) 
Eq.  (2.39) 
Eq.  (4.17) 
Eq.  (3.1) 

Eq.  (3.29) 
Eq.  (3.20) 
Eq.  (4.21) 
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Eq.  (4.18) 
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Eq.  (2.4) 
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Eq.  (3.17) 
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Eq.  (3.12) 

Eq.  (2.4) 
Eq.  (2.33) 
Eq.  (1.1) 
Eq.  (3.27) 


Fig.  3.4 
Fig.  3.3 


Fig.  3.3 
Fig.  2.1 


Fig.  3.4 
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T 

system  noise  temperature 

Eq.  (4.18) 

t 

time 

Eq.  (2.1) 

t 

running  variable  along  array 

Eq.  (3.28) 

U 

defocuslng  parameter 

Eq.  (3.15) 

ux,  u2 

generalized  defocuslng  parameters 

Eq.  (3-11) 

u 

Integration  parameter 

Eq.  (2.31) 

V 

integration  parameter 

Eq.  (2.31) 

X 

expansion  parameter 

Eq.  (2.21) 

xk 

projected  distance  to  kth  element 

Eq.  (3.34) 

Fig.  3.4 

y 

wavelength  directivity  parameter 

Eq.  (2.4?) 

yk 

projected  distance  to  kth  element 

Eq.  (3.35) 

Fig.  3.4 

yn 

imaginary  part  of  spherical  Bessel 
function  of  third  kind 

Eq.  (2.6) 

Greek 

a 

incidence  angle  of  line  of  sight 
from  transmitter 

Eq.  (1.1) 

Fig.  1.1 

3 

general  scattering  angle 

Eq.  (2.19) 

ei 

scattering  angle  of  first-order 
grating  lobe 

Eq.  (1.1) 

Fig.  1.1 

63dB 

3  dB  beamwith  in  angle 

Eq.  (2.49) 

r 

central  angle  between  users 

Eq.  (4.10) 

Y 

bistatic  scattering  angle  (also  6) 

Eq.  (4.9) 

Yk 

polar  angle  to  kth  element 

Eq.  (3.37) 

Fig.  3.3 

ab 

angular  departure  from  center  of 
first-order  grating  lobe 

Eq.  (2.43) 

A03dB 

3  dB  angular  departure  from  center 
of  first-order  grating  lobe 

Eq.  (2.49) 

AX 

wavelength  departure  from  center 
of  first-order  grat'Ug  lobe 

Eq.  (2.31) 

AA3dB 

3  dB  wavelength  departure  from 
center  of  first-order  grating 
lobe 

Eq.  (2.56) 

Af 

frequency  departure  from  center  of 
first-order  grating  lobe 

Eq.  (2.57) 

Af3dB 

3  dB  frequency  departure  from 
center  of  first-order  grating 
lobe 

Eq.  (2.58) 
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6k 

supplementary  scattering  angle  to 
kth  element 

Eq.  (3.20) 

e 

nadir  angle  of  line  of  sight 

Eq.  (3.19) 

C 

angle  between  line  of  sight  and 
plane  of  flexure 

Eq.  (3.43) 

n 

antenna  aperture  efficiency 

Eq.  (4.24) 

0P 

longitude  of  earth  point  (user) 

Eq.  (4.4) 

0s 

longitude  of  satellite 

Eq.  (4.4) 

• 

0s 

relative  orbital  angular  rate 
of  satellite 

Eq.  (4.3) 

0 

bistatic  scattering  angle 
(also  y ) 

Eq.  (2.3) 

0 

polar  angle  to  earth  reference 
point 

Eq.  (3.18) 

0 

max 

maximum  polar  angle  (to  edge  of 
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I.  SYNOPSIS 


The  survivability  and  endurance  of  the  communication  systems  used 
for  the  command  and  control  of  strategic  forces  are  vital  issues. 
Concern  that  existing  systems  may  be  vulnerable,  despite  their  diver¬ 
sity  and  specialized  design,  has  led  to  numerous  suggestions  for  their 
improvement.  Among  these  is  a  proposal  by  Yater  [1]  for  an  advanced 
passive  communication  satellite.  If  it  proves  practicable,  this 
design,  now  called  PACSAT,  may  be  useful  for  the  command  and  control 
of  strategic  forces  because  it  can  provide  survivable  point-to-point 
communication  at  low  to  moderate  data  rates  between  moderate  sized 
mobile  terminals  within  a  few  thousand  miles  of  one  another. 

A  typical  PACSAT  communication  link  is  illustrated  in  Fig.  1.1. 
PACSAT  consists  of  a  uniform  linear  array  of  scatterers  whose  axis  is 
pointed  (ideally)  toward  the  center  of  the  earth;  gravity-gradient 
stabilization  is  used  to  keep  the  array  erect  When  illuminated  by  a 
transmitter  on  the  earth,  the  array  scatters  energy  much  like  a  dif¬ 
fraction  grating.  The  element  spacing  is  chosen  such  that  only  the 
first  two  grating  lobes  are  generated.  Because  of  the  length  and  the 
radial  symmetry  of  the  array,  these  grating  lobes  are  in  the  form  of 
thin-walled  cones.  The  zeroth-order  grating  lobe  is  directed  out  to 
space  at  a  half  angle  a  equal  to  that  between  the  line  of  sight  from 
the  transmitter  to  the  array  and  the  array  axis.  The  first-order 
grating  lobe  is  directed  back  toward  the  earth  at  a  half  angle  B^, 
which  is  governed  by  the  relationship 

s(cos  a  +  cos  S^)  =  X  (1.1) 

where  s  is  the  element  spacing  and  X  is  the  operating  wavelength. 
Inasmuch  as  X  or,  equivalently,  the  frequency  f  can  be  chosen  to 
give  3^  any  desired  value,  the  scattered  beam  is  said  to  be  frequency- 
steerable. 

To  establish  a  communication  link  using  PACSAT,  an  operating  wave¬ 
length  is  chosen  such  that  the  annulus  illuminated  on  the  earth  by  the 
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the  two  lines  of  sight.  When  discussed  in  connection  with  scattering 
theory  this  angle  will  be  denoted  by  6. 

The  basic  command  and  contro  communication  problem  and  a  general 
description  of  the  PACSAT  concept  are  presented  in  a  companion  Note 
[2],  which  serves  as  a  comprehensive  introduction  to  this  report.  The 
companion  Note  can  profitably  be  read  first  to  establish  the  context 
for  PACSAT  and  to  gain  a  simple  understanding  of  its  basic  technical 
properties  and  its  system  characteristics  in  a  practical  application. 
Here,  the  electromagnetic  properties  and  communication  characteristics 
on  which  that  Note  draws  are  presented  in  detail.  The  present  report 
is  purely  technical  and  is  intended  only  to  provide  basic  design  in¬ 
formation  and  general  system  characteristics.  Issues  such  as  cost, 
jam  resistance,  physical  survivability,  erection  techniques,  replen¬ 
ishment  and  augmentation  strategies,  and  flexural  dynamics  will  be 
treated  in  forthcoming  publications. 

The  electromagnetic  properties  of  the  PACSAT  array  in  the  far 
field  are  treated  in  Sec.  II.  The  known  Mie  solution  [3]  for  scatter¬ 
ing  from  a  sphere  is  reviewed  and  adapted  for  the  case  of  circularly 
polarized  illumination  and  scattering,  which  is  the  one  of  practical 
interest.  The  gain  of  the  array  in  the  far  field  is  then  calculated 
assuming  uniform  illumination  along  the  array  (that  is,  end  effects 
are  disregarded).  These  results  are  combined  to  arrive  at  the  far- 
field  scattering  cross  section  of  the  array,  assuming  that  mutual 
effects  between  scatterers  can  be  ignored.  The  section  concludes  with 
a  computation  of  the  directivity  pattern  of  the  array  in  the  far 
field. 

Defocusing  caused  by  array  flexure  and  near-field  operation  are 
examined  in  Sec.  III.  It  is  assumed  that  the  effect  can  be  repre¬ 
sented  adequately  by  the  leading  nonlinear  term  in  the  expansion  of 
the  general  path  length  expression.  This  quadratic  term  is  then  used 
to  define  a  defocusing  parameter  that  completely  characterizes  the 
effect  of  array  flexure  and  near-field  operation  on  the  array  scatter¬ 
ing  cross  section  and  directivity  pattern.  The  specific  geometry  of 
PACSAT  in  orbit  is  then  used  to  relate  the  defocusing  parameter  to 
practical  array  and  system  parameters.  Results  are  presented  in  the 
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form  of  a  tradeoff  between  allowable  array  flexure  and  length  for  an 
acceptable  amount  of  defocusing  as  a  function  of  orbital  altitude  and 
array  libration  for  typical  user  locations.  Adequate  scattering  cross 
sections  (0.5  to  1.0  x  10^  m^)  can  be  realized  if:  (1)  the  orbital 
altitude  is  no  less  than  8000  km  or  so;  (2)  the  libration  is  limited 
to  2  or  3  deg  at  most;  and  (3)  the  flexure  does  not  exceed  an  eighth 
of  a  wavelength  or  so.  The  last  condition  is  particularly  stringent 
because  a  typical  operating  wavelength  may  be  only  3.75  cm  (8  GHz). 

The  system  implications  of  using  PACSAT  in  a  circular  equatorial 
orbit  are  considered  in  Sec.  IV.  As  noted  in  the  companion  Note,  this 
orbit  is  particularly  des_rable  because  it  not  only  limits  the  re¬ 
quired  beam-steering  frequency  excursions  and  facilitates  satellite 
acquisition  and  tracking,  but  also  provides  adequate  CONUS  coverage. 

The  geometrical  and  orbital  properties  are  investigated  to  show  how 
the  orbital  period,  relative  angular  rate,  and  coverage  are  related  to 
orbital  altitude.  Satellite  azimuth  and  elevation  angle  trajectories, 
orbital  arcs  of  mutual  visibility,  and  mutual  viewing  times  are  calcu¬ 
lated  for  typical  user  locations.  The  second  and  third  of  these  show 
that  satellite  spacings  of  about  60  deg  and  mutual  viewing  times  of 
about  60  min  can  be  achieved  only  by  the  use  of  orbital  altitudes  of 
8000  km  or  so.  Next,  some  major  array  operating  parameters  are  calcu¬ 
lated.  These  include  the  bistatic  scattering  angle  and  the  extremes 
of  operating  wavelength  that  would  be  experienced  by  typical  users  as 
a  function  of  orbital  altitude.  The  first  of  these  is  needed  to  cal¬ 
culate  the  scattering  cross  section  and  the  second  to  set  the  libra¬ 
tion  allowable  if  the  frequency  excursions  are  to  be  accommodated 
within  existing  frequency  allocations.  The  latter  leads,  typically, 
to  a  libration  limit  of  a  few  degrees  at  most. 

Representative  optimal  array  design  parameters  are  then  calculated 
and  it  is  seen  that  the  achievable  scattering  cross  section  ranges 
from  about  0.55L^  at  orbital  altitudes  of  4000  to  5000  km  to  about 
0.65L  at  10,000  km  (these  results  are  given  in  terms  of  the  array  length 
L).  Finally,  some  typical  communication  characteristics  are  presented 
in  the  form  of  power  budgets.  These  indicate  that  PACSAT  may  be  use¬ 
ful  for  low  data  rate  communication  between  mobile  or  transportable 


terminals.  For  example,  a  1  kbps  data  rate  can  be  supported  between 
1  m  diameter  dishes  at  a  frequency  of  8  GHz  using  a  10  kW  transmitter, 
a  low-noise  receiver,  and  a  PACSAT  array  with  a  scattering  cross  section 
of  10^  m^  (about  1300  m  long)  at  an  orbital  altitude  of  6000  km. 
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II.  ELECTROMAGNETIC  PROPERTIES  IN  THE  FAR  FIELD 


The  field  scattered  back  to  the  earth  by  the  PACSAT  array  at  a 
great  distance  is  investigated  in  this  section.  It  is  assumed  that 
PACSAT  consists  of  a  uniform  linear  array  of  perfectly  conducting 
spheres  whose  axis  is  pointed  (ideally)  toward  the  center  of  the 
earth.  First,  the  known  far-field  scattering  properties  of  a  single 
sphere  are  reviewed  and  its  scattering  cross  section  for  circularly 
polarized  illumination  and  scattering  is  calculated  over  the  range  of 
bistatic  scattering  angles  and  sphere  sizes  of  interest.  Then,  the 
gain  of  the  array  in  the  far  field  is  derived  and  used  to  calculate 
the  bistatic  far-field  scattering  cross  section  of  the  array. 

Finally,  its  directivity  pattern  and  beam  characteristics  in  the  far 
field  are  developed. 


FAR-FIELD  SCATTERING  FROM  A  SPHERE 

The  problem  of  scattering  from  a  sphere  was  solved  exactly  by  Mie 
[3]  in  1908  and  Is  discussed  in  Stratton  [4];  the  results  presented 
here  are  from  Barrick  [5J. 

Consider  an  Incident  plan  wave 


-i  i  'i(koz+V>- 

E  =  E  e  x 

x  x 


(2.1) 


where  x  is  a  unit  vector  in  the  x  direction. 


k 


0 


(2.2) 


and  X  is  the  wavelength;  the  harmonic  time  dependence  will  be  sup¬ 
pressed  henceforth.  As  illustrated  in  Fig.  2.1,  the  incident  field  is 
linearly  polarized  in  the  x-z  or  E-plane  and  is  propagating  in  the  -z 
direction.  A  sphere  of  radius  a  is  at  the  origin  and  the  scattered 
field  is  desired  at  a  point  at  a  distance  r,  polar  angle  0,  and 
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Fig.  2.1  —  Geometry  of  field  scattered  from  a  sphere 


meridional  angle  <f> ,  with  respect  to  the  x-z  plane.  The  angle  0  is 
also  known  as  the  bistatic  scattering  angle;  in  the  monostatic  case, 
0=0  and  the  field  is  said  to  be  backscattered . 


*The  bistatic  scattering  angle  is  denoted  by  0  in  this  section  to 
conform  with  standard  usage  in  scattering  theory.  Inasmuch  as  0  is 
also  conventionally  used  to  denote  the  polar  angle  in  spherical  geome¬ 
try,  it  is  used  for  that  purpose  elsewhere  in  this  report  and  y  is 
then  used  to  denote  the  bistatic  scattering  angle. 
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When  the  distance  is  great,  the  scattered  field  is  essentially  a 
spherical  wave  that  can  be  expressed  as 


ikQr 

E®(6  ,<*>)  =  S__  |Sl(6)  cos  4>  6  -  S2(0)  sin 


(2.3) 


The  functions  and  S2  are  given  by 


s1(©>  =  X  (-D 


n+1 


P  (cos  0 )  a  l 

A  — 5 +  i  B  ^-P  (cos  0) 

n  sin  6  n  do  n 


(2.4) 


oo 

s,(0)  =  X)  (-1) 


n+1 


n=l 


P*(cos  0) 

A  ~  P  (cos  6)  +  i  B  -2— — s— 
n  d6  n'  n  sin  0 


where  P^  is  the  associated  Legendre  function.  For  a  perfectly  con¬ 
ducting  sphere  of  radius  a,  the  coefficients  An  and  Bn  are  given  by 


.  nn  2n  +  1  jn(k0a) 
n  "  n(n  +  1)  .  (1) 


hn  (k0a) 


(2.3) 


o  f  nn41  2n  +  1  lkOaVkOa> 

"  '  +  U  >V  h'l>(koa) 


where  the  prime  denotes  differentiation  with  respect  to  kQa  and 

hn1)(z)  “  jn(z)  +  1  yn(z)  (2'6) 


is  the  spherical  Bessel  function  of  the  third  kind. 
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The  8  and  $  components  of  the  scattered  field  are  depicted  in 
Fig.  2.1.  It  is  apparent  from  Eq.  (2.3)  that  there  is  only  a  9  compo¬ 
nent  in  the  E-plane  (<t>  =  0)  and  only  a  $  component  in  the  H-plane 
(<{>  =  tt / 2 ) ;  thus,  the  scattered  field  is  linearly  polarized  in  these 
planes.  Elsewhere,  because  these  components  are  complex  and  unequal, 
in  general,  the  scattered  field  is  elliptically  polarized  (though  only 
weakly  so  for  the  range  of  values  of  kQa  and  8  of  interest). 

A  more  general  representation  in  which  the  incident  linearly 
polarized  field  has  an  arbitrary  orientation  can  be  obtained  by  con¬ 
sidering  a  second  incident  field  that  is  linearly  polarized  in  the  y 
direction.  Thus,  for 


pi  "lkoz- 

E  e  y 

y 


(2.7) 


the  scattered  field  can  be  written  by  substituting  ij>  -  tt/2  for  in 
Eq.  (2.3)  to  obtain 

]  (2.8) 


ikQr 

El(9’*)  =  Ey  V 


S L (0  )  sin  ij>  8  +  (Q  )  cos 


The  arbitrarily  oriented,  linearly  polarized,  incident  plane  wave 
is  then  equal  to  the  sum  of  Eqs.  (2.1)  and  (2.7),  or 


E1  =  E1  +  E1 
x  y 


(2.9) 


and  the  general  total  scattered  field  is  given  by  the  sum  of  Eqs. 
(2.3)  and  (2.8),  or 
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Es(e  ,<t>) 


eJ(0,$)  e  +  eJ(0 ,♦>  ♦ 

^  [  /i  i  v. 

Vt"  lSl(e)(Ex  cos  ♦  +  eJ  sin 

+  S2(0  )  E*  sin  4>  +  E*  cos  <J>)  ? 


(2.10) 


In  matrix  form, 


where 


this  becomes 

Eq(6,*) 

E|<8  ,<}>) 


rn 

a6x  a0  y 

E 

X 

1 

a,  a. 

E 

<j>x  <(>y 

y 

(2.11) 


ik0r 


ik0r 


6x  kQr  “1 


ik0r 


S  (6)  cos  $ , 


6 y  kQr  -1 


S.  (0  )  sin  <J> 


(2.12) 


ik0r 


a<(»x  =  -~k7-S2(0)  8ln  *’ 


<t>y  ^0r  2 


S_(0  )  cos  <f> 


The  double  subscript  notation  on  the  coefficients  of  the  scattering 
matrix  refers  to  the  polarization  of  the  scattered  and  incident  field 
components,  in  that  order. 

In  the  PACSAT  application,  it  is  desirable  to  use  circular  rather 
than  linear  polarization  for  both  transmission  and  reception  because 
it  eliminates  the  need  for  polarization  alignment.  The  representation 
of  the  scattered  field  in  terms  of  right-  and  left-hand  circularly 
polarized  components  when  the  incident  field  has  both  right-  and  left- 
hand  circularly  polarized  components  can  be  written  in  matrix  form  in 
a  fashion  analogous  to  Eq.  (2.11),  or 
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*J(M> 

aRR 

aRL 

Ej(e^) 

s 

(2.13) 

aLR 

aLL_ 

LelJ 

The  relationships  between  the  elements  of  the  circular  and  linear 
polarization  matrices  are  given  by  Barrick  (Ref.  5,  Eqs.  (2.1-25))  as 


aRR  -V2((-  %y  +  a0x)~  1  (a8y  +  v)|-  aRL  “  ^|f%,  +  a0x)+  1  (*8 y  '  %xj) 

(2 . 14) 

aLR  “^|Ky  +  a9x)  *  1  fa6y  '  %*))•  aLL  "  a»y  +  a8x)  +  1  fa9 y  +  V)] 


Substituting  Eqs.  (2.12)  in  Eqs.  (2.14)  then  yields 


RR 


ik0r 


s^e)  -  s2(9)|  cos 


i 


s1(6)  -  S2(0) 


(2.15) 


LR 


~2k~r  |{Si(0  )  +  S2(0)]cos  *  -  i  ^(0)  +  S2(0  )]  sin  $J 


It  can  be  seen  from  Eqs.  (2.14)  that  the  other  coefficients  are  given 
by 


LL 


aRR  6 


12<|> 


aRL°  aLR® 


i2<J> 


(2.16) 


Finally,  the  various  scattering  cross  sections  are  given  by 
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o 


RR 


a 


RL 


(2.17) 


so  substituting  from  Eqs.  (2.15)  and  (2.16)  yields 


a  =o 

RR  LL 


-§-  |  [Re  S^e)  -  Re  S2(0  )j 2  +  [  Im  S^B)  -  1m  S2(6  )J2j 


(2.18) 


°LR  aRL  ,2 
k0 


Re  S.(0)  +  Re  So(0)f  +  Im  S.  (6  )  +  Im  S„(6) 


where  Re  and  Im  denote  the  real  and  imaginary  parts,  respectively. 

For  the  PACSAT  application,  the  bistatic  scattering  angle  6  will 
generally  be  less  than  about  45  deg  or  so.  Furthermore,  to  maximize 
the  field  scattered  toward  the  receiver,  it  is  desirable  to  operate  in 
the  vicinity  of  the  first  resonance,  which  occurs  at  values  of  k0a  of 
about  1.0  to  1.1.  (It  will  be  seen  later  that  the  peak  of  the  reso¬ 
nance  for  the  array  of  spheres  shifts  the  values  of  kga  to  about  1.1 
to  1.2.)  Therefore,  the  scattering  cross  sections  have  been  evaluated 
numerically  for  kga  values  of  between  0.9  and  1.4  and  values  of  0  from 
0  to  50  deg.  The  same-sense  scattering  cross  section  o  or  o  is 

Kiv  LL 

small  and  of  no  practical  interest.  The  opposite-sense  scattering 

cross  section  a  or  a  ,  which  conveys  virtually  all  of  the  scattered 
LK  KL 

energy,  is  plotted  in  Fig.  2.2  after  normalization  to  its  physical 
2 

cross  section  na  .  (See  the  appendix  for  details  of  the  computation.) 

To  facilitate  reading  the  three-dimensional  plot,  contours  of 
2 

constant  oLR/na  are  shown  as  solid  lines  on  the  surface.  Also  shown, 
as  dashed  lines,  are  contours  of  constant  0  and  kga.  A  heavy  dashed 
line  indicates  the  ridge  along  which  the  function  takes  on  a  local 
maximum.  As  mentioned  earlier,  the  contour  for  8  *  0  is  the 
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Normalized 
backscattering 
cross  section 
a 

na1 


Fig.  2.3  —  The  normalized  radar  or  backscattering  cross  section 
of  a  perfectly  conducting  sphere 


GAIN  OF  THE  ARRAY  IN  THE  FAR  FIELD 

The  far-field  geometry  of  the  PACSAT  array  is  shown  In  Fig.  2.4. 
The  N  scattering  elements,  which  are  assumed  to  be  isotropic,  are 
spaced  a  uniform  distance  s  along  the  array  axis.  When  the  array  is 
illuminated  from  a  great  distance  at  an  angle  a  with  respect  to  its 
axis,  the  incoming  rays  can  be  considered  parallel.  The  rays  repre¬ 
senting  the  scattering  toward  a  distant  point  at  an  angle  (5  can  then 
also  be  considered  parallel.  It  is  clear,  from  symmetry,  that  the 
planes  containing  the  incident  and  scattered  rays  need  not  be  coinci¬ 
dent,  i.e.,  there  is  axial  symmetry. 

The  uniform  element  spacing  results  in  a  uniform  path  lengthening 
or  phase  progression  between  successive  rays.  It  will  be  assumed  that 
the  elements  are  uniformly  illuminated,  which  is  equivalent  to  ignor¬ 
ing  end  effects.  It  is  well  known  that  the  currents  induced  on  pas¬ 
sive  array  elements  are  influenced  by  the  presence  of  other  nearby 
elements.  When  the  array  symmetry  is  appropriate,  as  when  the  elements 
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N  scattering 


axis 


Fig.  2.4  —  The  basic  PACSAT  array 


are  In  a  ring,  these  mutual  effects  will  be  the  same  for  all  elements 
and  there  is  no  approximation.  If  there  are  a  large  number  of  ele¬ 
ments  in  non-symraetrical  configurations,  as  would  be  the  case  for 
PACSAT,  only  the  elements  near  the  ends  of  the  array  would  experience 
significantly  different  mutual  effects.  This  end  effect  is  considered 
negligible  for  PACSAT  and  will  be  ignored. 

The  magnitude  of  the  scattered  field  can  be  written 


. N-l  I 

|e(  -  I  £  exp  i  ns(cos  a  +  cos  8)J 

n-0 


(2.19) 
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where  the  Individual  elements  are  assumed  to  contribute  unit  amplitudes 
at  the  receiving  point.  The  sum  can  be  expressed  in  closed  form  as 


where 


sin  |nx 

~  I 

sin  2* 


(2.20) 


x  =  2n  y  (cos  a  +  cos  3)  (2.21) 

Inasmuch  as  the  spacing  s  and  the  illumination  angle  os  are  parameters 
of  the  problem,  the  amplitude  of  the  scattered  field  is  a  function  of 
two  variables:  the  scattering  angle  0  and  the  wavelength  X.  Also, 
the  scattered  field  has  maxima,  of  amplitude  N,  called  grating  lobes, 
whenever 


j-  (cos  a  +  cos  0)  =  k,  k  =  0,  1,  2,  ...  (2.22) 

where  k  is  the  order  of  the  grating  lobe. 

The  zeroth-order  grating  lobe  (k  =  0)  occurs  when 

0  =»  n  -  a  (2.23) 

so  it  is  independent  of  X ;  it  is  shown  in  Fig.  1.1  as  the  grating  lobe 
scattering  off  into  space  (for  practical  values  of  a).  The  first- 
order  grating  lobe  (k  =  1)  occurs  when  6  and  X  satisfy 

s(cos  a  cos  0)  -  X  (2.24) 

(Grating  lobes  of  higher  order  are  suppressed  in  the  PACSAT  applica¬ 
tion  by  avoiding  values  of  X  small  enough  to  satisfy  Eq.  (2.22)  for 
k  >  2.)  It  is  convenient  to  specify  the  center  of  the  first-order 
grating  lobe  given  by  Eq.  (2.24)  by  either  fixing  X  and  letting 
8  •  8^  denote  the  angle  at  which  it  occurs,  i.e., 

s(cos  a  +  cos  6^)  -  X 


(2.25) 
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or  fixing  8  and  letting  X  =  X^  denote  the  wavelength  at  which  It 
occurs,  i.e., 


s(cos  a  +  cos  8)  =  X^  (2.26) 

The  first  of  these  is  used  in  Eq.  (1.1)  and  in  Fig.  1.1  to  depict  the 
frequency-steerable  first-order  grating  used  to  form  a  communication 
link.  It  will  also  be  used  in  this  subsection  to  integrate  the  scat¬ 
tered  field  over  8.  Both  will  be  used  later  to  describe  the  directiv¬ 
ity  pattern  in  the  vicinity  of  the  main  lobe. 

The  gain  of  the  array  is  defined  as  the  ratio  of  the  field  inten¬ 
sity  in  the  direction  of  maximum  radiation  to  the  average  intensity  of 
radiation  in  all  directions.  Thus, 


E 

2 

max 

• 

E 

2d0 

(2.27) 


where  the  integral  is  over  the  sphere.  It  is  apparent  from  Eq.  (2.20) 
that 


Ul 


2 

1  max 


(2.28) 


To  facilitate  the  integration  on  8,  use  Eq.  (2.25)  to  eliminate  X  in 
Eq.  (2.21).  Then,  using 


2tt  c03  a  +  cos  3 
cos  a  +  cos  8^ 


(2.29) 


as  a  change  of  variables  in  Eq.  (2.20),  and  conventional  spherical 
variables,  the  integral  in  the  denominator  of  Eq.  (2.27)  becomes 
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„  ,  2  1„ 

2n  it  sin  -Nx 

J  =  f  d<j>  C  d8  - r — -  sin  8 

JO  Jo  sin  |x 

,  2  K 

u  sin  xNx 

=  |cos  a  +  cos  $])f  - 2~ I —  dx  (2.30) 

A  sin  2X 


where 


u  =  2it 


cos  a  +  1 
cos  a  +  cos  8, 


2r 


cos  a  -  1 
cos  a  +  cos  8 


(2.31) 


1 


The  integrand  of  Eq.  (2.30)  is  plotted  in  Fig.  2.5  for  N  =  12  and 

17  and  a  -  8^  =  30  deg  to  illustrate  its  behavior  (the  abscissa  is 

2 

linear  in  x).  The  function  attains  its  maximum  value  of  N  at  x  =  0, 
2n,  which  correspond  to  the  values  of  8  for  the  zeroth-  and  first- 
order  grating  lobes,  respectively.  Even  for  small  numbers  of  elements 
in  the  array,  the  principal  contribution  to  the  integral  is  seen  to  be 
in  the  vicinity  of  these  values  of  x.  As  N  increases,  the  peak  value 
increases  rapidly  and  the  integrand  goes  through  an  increasing  number 
of  oscillations  (these  are  actually  the  sidelobes  of  the  array  pattern 
in  the  angle  8).  Also,  the  integrand  compacts  about  the  values  at 
which  it  attains  its  peak  value.  Inasmuch  as  practical  values  of  N 
will  be  in  the  tens  of  thousands,  it  is  clear  that  the  integral  can  be 
well  approximated  by  writing  it  as  the  sura  of  two  overlapping  infinite 
integrals  in  which  the  sine  function  in  the  denominator  is  expanded 
about  x  *  0,  2x .  Thus  Eq.  (2.30)  can  be  written 


2  1  2  1 
“  sin  rNx  u  sin 

J  -  cos  a  +  cos  8.  f  — j ~ —  dx  +  f  -rj - =-  dx 

Jy  (fx>  •/—  l|(x  -  2n)r 

(2.32) 


The  limits  u  and  v  have  been  retained  In  Eq.  (2.32)  to  accommodate  the 
cases  when  v  is  near  0  and  u  is  near  2tt  ,  which  correspond  to  a  and 
8^  small,  respectively. 


V  0 

(=  -0.4860) 


2n  u 
(=  6  7692) 

—4 - 1 

6,  0 


p  - — I — I - 

n  n-a 

(a)  1 2  elements  (N  =  12) 


Fig.  2.5  —  Plots  of  integrand  of  Eq.  (2.32)  for  a  =0,  =  n/ 6 


The  integrals  in  Eq.  (2.32)  are  easily  evaluated  and,  when  sub¬ 
stituted  in  Eq.  (2.27),  yield 


-20- 


2nN 


(cos  a  +  cos 


V 


+  si  N(u  -  2n) 


1  -  cos  N(u  -  2n  ) 
N(u  -  2n  ) 


-  8i  Nv  + 


1  -  cos  Nv 
Nv 


(2.33) 


where  Si  is  the  sine  integral.  Limiting  cases  occur  when:  (1)  both 
a  and  0^  are  equal  to  zero;  (2)  when  either  a  or  0  ^  is  equal  to  zero 
and  the  other  is  sufficiently  large  that  its  associated  integral  in 
Eq.  (2.32)  takes  on  its  value  for  infinite  limits;  and  (3)  when  both 
a  and  0^  are  similarly  sufficiently  large.  To  establish  such  a  cri¬ 
terion  for  a  and  0^,  recall  that  the  principal  contribution  to  the 
integrals  in  Eq.  (2.32)  comes  from  the  main  lobe  of  the  integrand.  It 
can  be  shown  that  if  the  integration  is  to  the  first  null,  the  inte¬ 
gral  attains  95  percent  of  its  final  value;  if  the  integration  is  to 
the  fifth  null,  the  integral  reaches  99  percent  of  its  final  value. 
Thus,  if  a  (or  0^)  is  not  "small,"  where  "small"  means  on  the  order  of 
a  few  beamwidths,  its  associated  integral  in  Eq.  (2.32)  essentially 
takes  on  the  value  for  v  =  (or  u  =*  -t*° )  •  The  array  gain,  Eq. 

(2.33),  can  then  be  specialized  to  the  following  cases: 


G 


N, 


1  +  cos  0^’ 


1  +  cos  a’ 
N 


cos  a  +  cos  0^’ 


a  “  0  »■  0 

a  *  0,  0^  not  small 

0^  =  0,  a  not  small 
a ,  0  not  small 


(2 .34) 


The  cases  for  either  a,  0  ,  or  both  small  or  equal  to  zero,  all  of 
which  result  in  end-fire  operation  of  the  array,  are  not  of  practical 
interest  because  they  correspond  to  fleeting  or  unlikely  geometries. 
Most  of  the  time,  the  gain  will  be  given  by  the  fourth  of  Eqs. 

(2.34).  Inasmuch  as  it  is  also  a  conservative  value,  it  will  be  used 
in  the  following  analysis. 
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FAR-FIELD  SCATTERING  FROM  THE  ARRAY 

The  scattering  cross  section  a  of  the  array  Is  defined  as 

o  =  GA  ,,  (2.35) 

ef  1 

where  G  Is  the  array  gain  given  by  Eqs.  (2.34)  and  Ae££  is  its  effec¬ 
tive  capture  cross  section.  As  noted  above,  the  gain  will  be  taken  as 
the  fourth  of  Eqs.  (2.34).  It  will  be  assumed  that  mutual  effects 
between  the  Individual  spheres  comprised  by  the  array  can  be  ignored 
so  that  the  effective  capture  cross  section  of  the  array  can  be  taken 
as  N  times  the  scattering  cross  section  of  an  isolated  sphere,  i.e., 


A 


eff 


No 


LR 


(2.36) 


where  a  is  given  by  Eq.  (2.18).  The  magnitude  of  the  error  inherent 
LR 

in  the  assumption  is  not  known  but,  unfortunately,  an  exact  analysis 
does  not  appear  possible  so  there  is  no  choice.  It  is  hoped,  by  anal¬ 
ogy  to  the  known  properties  of  conventional  antenna  arrays,  that  the 
error  is  not  significant. 

Substituting  from  the  fourth  of  Eqs.  (2.34)  and  using  Eq.  (2.25) 
to  reintroduce  X  leads  to 


a 


N2so 


LR 


(2.37) 


I'^ing  Eq.  (2.2),  to  introduce  kg,  this  can  be  brought  into  the  form 


0  -^-f)koa^  (2-38> 

2 

which  permits  identifying  N  sa/2  -.s  an  appropriate  normalization  fac- 

2 

tor.  The  term  a  /  a  is  the  normalized  scattering  cross  section  of  a 
LR 

sphere  plotted  on  Fig.  2.2,  so  t’..i  numerical  values  need  only  be 
multiplied  by  kga  to  obtain  the  values  of  the  normalized  array 
scattering  cross  section. 
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Other  forms  of  the  normalization  factor  can  be  obtained  by  noting 
that  the  array  length  L  Is  given  by 


L  s  (N  -  l)s  2  Ns 


(2.39) 


Also,  from  Eq.  (2.25),  an  end-fire  operating  wavelength  AQ  can  be 
identified  by  setting  a  =  0^  =  0  to  get 


Fig.  2.6  —  Scattering  from  the  PACSAT  array  of  length  L  at  wavelength  A 


With  the  help  of  these,  it  follows  that 


N2sa 

2 


N2aX , 


(2.41) 


These  are  equivalent  and  may  be  interchanged  as  preference  dictates. 

Using  the  simplest  form,  the  normalized  array  scattering  cross  section 
2 

oAq/L  a  is  plotted  in  Fig.  2.6  in  a  form  similar  to  that  used  in  Fig. 
2.2.  The  value  of  k.Qa  at  which  the  maximum  is  attained  is  plotted  in 
Fig.  2.7  to  facilitate  selecting  the  sphere  size  that  optimizes  the 
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Fig.  2.7  —  Optimum  scattering  from  the  PACSAT  array 
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design  of  the  array.  A  practical  design  procedure  will  be  outlined  in 
the  next  section  and  an  example  shown  for  PACSAT  in  a  circular  equa¬ 
torial  orbit. 


DIRECTIVITY  PATTERN  IN  THE  FAR  FIELD 

The  far-field  directivity  pattern  D  is  given  by  the  ratio  of  the 
square  of  Eq.  (2.20)  to  Eq.  (2.28)  so 

/sin  |nx  V 

D  =  - -  (2.42) 

\N  sin 

where  x  is  given  by  Eq.  (2.21) 

x  =  2ir  (cos  a  +  cos  0)  (2.21) 

The  directivity  pattern  is  desired  in  the  vicinity  of  the  first-order 
grating  lobe.  Inasmuch  as  this  grating  lobe  is  a  function  of  both 
angle  and  frequency,  it  is  most  conveniently  expressed  by  two  compo¬ 
nent  directivity  patterns.  One  is  the  conventional  pattern  in  which 
the  wavelength  is  fixed  and  the  pattern  is  expressed  in  terms  of  the 
angular  departure  from  the  angle  corresponding  to  the  peak  of  the 
first-order  grating  lobe.  In  the  other,  the  angle  is  fixed  and  the 
pattern  is  expressed  in  terms  of  the  wavelength  departure  from  the 
wavelength  corresponding  to  the  peak  of  the  first-order  grating  lobe. 


Angle  Pattern 

For  this  pattern  it  is  appropriate  to  use  Eq.  (2.25)  to  eliminate 
X  in  Eq.  (2.21),  so,  as  above, 


cos  a  +  cos  8 
cos  a  +  cos  6  ^ 


(2.29) 


Let  the  angle  B 
beam  center  B^, 


be  expressed  in  terms  of  its  departure  A0  from  the 
i .  e . , 
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@  -  01  +  A0 


(2.43) 


Then,  for  Ag  small 


cos  0  *  cos  0^  -  A0  sin  0^ 


(2.44) 


and 


sin  =  ±  sin 


NttA0  sin  0 . 


I  cos  at  +  cos  0, 


(2.45) 


where  the  signum  is  -  for  N  even  and  +  for  N  odd.  For  N  =  1  and  A0 
small , 


nA0  sin  0, 


sin  rx  = 


2  cos  a  +  cos  0, 


(2.46) 


so  substituting  Eqs.  (2.45)  and  (2.46)  in  Eq.  (2.42)  yields 


”  -  M 


(2.47) 


where 


NttA0  sin  0 


1 


cos  a  +  cos  0 


1 


(2.48) 


The  3  dB,  or  half-power,  point  occurs  when  y  =  1,39156,  so  the  3  dB 
beamwidth  is  given  by 


cos  a  +  cos  0 

0,.n  =  2A0_  -  0.8859  - rr — t — s - - 

3dB  3dB  N  sin  0 


(2.49) 
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s i.n  y 

The  result  is  invalid  for  0  near  zero  where  the  - *-  patterns  char- 

^  s  i  n  y 

acterizing  the  walls  of  the  conical  pattern  merge  into  a  single  — — *- 
pattern  characterizing  the  end-fire  pencil  beam. 

The  near-inverse  relationship  between  the  3  dB  beamwidth  0^^ 
and  the  cone  half  angle  0^  in  Eq.  (2.49)  indicates  how  the  wall  thick¬ 
ness  of  the  scattered  cone  decreases  as  the  cone  opens.  Note  that  the 
total  solid  angle  in  the  cone  is  approximately  equal  to  2n0.jdB  sin  0  ^ . 
Thus,  the  total  solid  angle  is  roughly  constant;  this  is  consistent 
with  Eqs.  (2.34),  which  show  that  the  gain  is  also  roughly  constant. 
The  normalized  3  dB  beamwidth,  is  plotted  as  a  function 

of  0j  in  Fig.  2.8  for  a  range  of  illumination  angles  a. 

Wavelength  Pattern 

For  this  pattern  it  is  appropriate  to  use  Eq.  (2.26)  to  introduce 
X^  in  Eq.  (2.21)  to  obtain 

X1 

x  =  2*  -y  (2.50) 

Let  the  wavelength  X^  be  expressed  in  terms  of  its  departure  AX  from 
the  beam-center  value  X^,  so 


X  =  +  AX 


Then,  for  AX  «  X^ 

which  leads  to 


sin  ^Nx 


±  sin  Nn 


AX 

X, 


where  the  signum  is  -  for  N  even  and  +  for  N  odd,  and 


sin  2* 


(2.51) 


(2.52) 


(2.53) 


(2.54) 


2  3  4  5  7  10  20  30  40  50  70 

Half  angle  of  scattered  cone  /3i.  deg 
Fig.  2.8  —  Normalized  3  dB  beamwidth  of  PACSAT  array 


Substituting  Eqs.  (2.53)  and  (2.54)  In  Eq.  (2.42)  then  yields 


100 


/sin  Nh 
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The  3  dB,  or  half-power,  point  occurs  when  the  argument  equals 
1.39156,  so  the  3  dB  beamwidth  is  given  by 


*X3dB  '  °‘8859  T 


(2.56) 


The  basic  pattern,  including  the  first  five  sidelobes,  is  plotted 
in  Fig.  2.9  as  a  function  of  angle,  Eq.  (2.47),  wavelength,  Eq.  (2.55), 
and  frequency,  Eq.  (2.57).  A  numerical  example  helps  to  give  a  feel 
for  the  beamwidths  involved.  Thus,  for  N  =*  80,000  and  a  -  3^  *  15 
deg,  ^3dB  *  8.266  x  10  ^  rad  =  0.00474  deg;  if  a  ■  8  ^  *  30  deg,  * 

3.836  x  10  ^  rad  =  0.00220  deg.  The  extreme  thinness  of  the  cone  of 
scattered  radiation  is  apparent;  the  practical  concomitant  is  an  un¬ 
precedented  precision  required  in  beam  pointing. 

Again,  for  N  *  80,000  and  f^  *  8  GHz,  f 3^^  =*  88.59  kHz.  This 
result  highlights  the  inherent  narrow-bandedness  of  the  array.  First, 
it  gives  a  measure  of  the  frequency  control  required  because  it  shows, 
in  this  case,  that  a  frequency  change  of  only  88.59  kHz,  which  is 
roughly  one  part  in  10^  of  the  radiated  frequency,  will  steer  the 
scattered  beam  from  one  3  dB  point  through  the  peak  of  the  beam  to  the 
other  3  dB  point.  Second,  it  indicates  that  modulation  should  be  such 
that  the  modulated  signal  spectrum  not  significantly  exceed  the  3  dB 
bandwidth  of  the  array  because  the  array  is  acting  as  a  band-pass 
filter  that  will  attenuate  the  distant  sidebands. 
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Normalized  pattern  variable 
NA>3  sin 

cos  a  +  cos/3,  N  A,  N  f. 


Fig.  2.9  —  Far-field  directivity  pattern  of  PACSAT  array 
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III.  QUADRATIC  PHASE  ERROR  EFFECTS 


If  PACSAT  were  perfectly  straight  and  at  a  great  distance,  the 
properties  described  in  the  preceding  section  would  apply.  However, 
these  conditions  do  not  obtain  in  practice  and  it  is  necessary  to 
estimate  the  effects  of  nonideal  operation.  In  this  section,  a  quad¬ 
ratic  phase  error,  which  is  assumed  to  be  characteristic  of  small 
departures  from  ideal  conditions,  is  introduced  to  display  its  effect 
on  the  directivity  pattern  and  the  scattering  cross  section.  These 
results  are  then  related  to  the  specific  effects  associated  with  near¬ 
field  operation  and  flexure  of  the  array. 

THE  PEFOCUSING  PARAMETER 

In  Eq.  (2.19)  the  amplitude  of  the  distant  field  was  written  in 
terms  of  the  phase  progression  associated  with  the  uniformly  increas¬ 
ing  path  lengths  via  the  various  scatterers.  To  accommodate  nearby 
terminals  or  a  nonstraight  array,  the  general  form 


where  is  the  path  length  from  the  transmitter  to  the  receiver  via 
the  kth  element,  is  needed.  The  path  length  can  be  expanded  in  a 
Taylor's  series 


(3.2) 


about  some  convenient  point,  which  will  be  taken  here  as  the  center  of 
the  array. 

The  initial  constant  Cq  is  readily  identified  as 


(3.3) 
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where  Lq  Is  the  path  length  from  the  transmitter  to  the  receiver  via 
the  center  element  of  the  array.  The  first  coefficient  c^  gives  the 
linear  progression  In  path  length  used  In  Eq.  (2.19).  Inasmuch  as 
Interest  here  Is  centered  on  the  behavior  of  the  first-order  grating 
lobe,  this  coefficient  will  be  taken  as 

Cj  *  \  (3.4) 

where  X^  Is  the  wavelength  at  which  the  peak  of  this  grating  lobe  Is 
formed. 

Through  terms  of  second  order, 

Lk  ”  L0  +  Xlk  +  c2k2  (3.5) 

Taking 

X  »  Xx  -  AX  (3.6) 

then  leads,  for  AX  «  X^,  to 


Substituting  Eq.  (3.7)  In  (3.1)  and  discarding  terms  of  unit  magnitude 
results  In 


N-l 

2 


exp  1  2* 


(3.8) 


where  the  number  of  elements  N  must  be  odd  to  place  one  at  the  center  , 

of  the  array. 

For  N  large,  the  sum  can  be  approximated  by  an  integral  f 


\ 


-32- 


N-l 


Ei  - 1 

f  exp  i  2tt  | 

AA  C2  2\ 

~  x  +1~  x 

J. 

A  1  / 

2 

Completing  the  square  and  Integrating  leads  to 


(3.9) 


(3.10) 


where 


U 


1 


(3.11) 


When  written  with  a  lower  limit  of  zero,  the  integral  in  Eq.  (3.10) 
has  real  and  imaginary  parts  given,  respectively,  by  the  cosine  and 
sine  Fresnel  Integrals 


C(z) 


2  Z 

J  cos  J  t2  dt,  S(z)  -  j  sin  ~  t2  dt  (3.12) 


Using  these  functions,  the  normalized  directivity  pattern  becomes 


N2  I 


|c(u1) 


C(U2) 


sCUj)  -  S(U2) 


(3.13) 
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At  beam  center  A  A  »  0,  and 


where,  for  N  large, 


(3.14) 


(3.15) 


Inasmuch  as  and  U£  are  easily  expressed  in  terms  of  U  by 

naa/a  naa/a 

U2'-U - +U  <3’16> 

it  is  seen  from  Eq.  (3.15)  that  U,  which  will  be  called  the  defocusing 
parameter,  completely  characterizes  the  effect  of  second-order  phase 
errors  on  the  array  performance.  This  is  illustrated  in  Fig.  3.1, 
which  shows  the  mainlobe  and  the  first  five  sidelobes  of  the  direc¬ 
tivity  pattern  given  by  Eq.  (3.13)  for  representative  values  of  U.  I 

The  far-field,  or  focused  pattern,  which  is  identical  to  the  one 
plotted  in  Fig.  2.9,  results  when  U  ■  0.  As  the  defocusing  parameter 
is  increased,  the  directivity  (and,  hence,  the  array  scattering  cross 
section)  falls  off  as  the  sidelobe  levels  are  raised  and  their  struc¬ 
ture  is  altered.  For  U  -  1.0,  the  scattering  cross  section  is  reduced 
by  about  1  dB,  the  first  null  is  virtually  eliminated,  and  the  first 
sidelobe  is  raised  by  about  3  dB.  The  effect  on  higher-order  side- 
lobes  is  not  significant  beyond  about  the  third  sidelobe. 

Further  increases  in  the  defocusing  parameter  depress  the  center 
of  the  main  lobe,  rapidly  at  first,  then  in  a  steadily  increasing 

oscillatory  fashion.  (This  is  shown  in  Fig.  3.2,  which  is  a  plot  of  , 

Eq.  (3.14).)  At  the  same  time,  the  first  few  sidelobes  are  radically 

, 
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Fig.  3.1  —  Directivity  pattern  of  defocused  PACSAT  array 


r 
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Fig.  3.2  —  Scattering  cross  section  loss  for  the 
defocused  PACSAT  array 

affected.  For  example,  when  U  -  2.0,  the  scattering  cross  section 
loss  Is  about  10  dB  and  the  first  three  or  four  sldelobes  have  merged 
Into  a  single  lobe  with  a  3  dB  dip  at  Its  center  and  a  3  dB  beamwidth 
about  five  times  as  great  as  that  for  the  far-field  pattern.  All 
semblance  to  a  normal  directivity  pattern  in  the  vicinity  of  the  main 
lobe  la  lost  by  the  time  U  ■  3.0.  However,  even  large  values  of  U 
have  little  effect  on  sldelobes  that  are  well  removed  from  the  main- 
lobe.  For  Instance,  when  U  ”  5.0,  the  100th  sidelobe  has  its  peak 
value  raised  by  only  about  0.5  dB;  the  main  effect  Is  to  fill  in  the 
nulls  somewhat.* 


^Although  the  100th  sidelobe  may  seem  distant  In  terms  of  ordi¬ 
nary  antenna  practice,  it  should  be  noted,  from  Eq.  (2.42),  that  the 
PACSAT  directivity  pattern  has  N  lobes. 
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The  general  properties  of  the  defocusing  parameter,  as  it  is 
related  to  the  system  parameters,  are  developed  and  discussed  below. 
Although  the  choice  is  necessarily  subjective,  a  value  of  U  *  1.0  will 
be  taken  here  as  the  desirable  maximum  value  of  the  defocusin-;  param¬ 
eter.  It  will  be  used  in  the  next  section  to  develop  a  tradeoff  be¬ 
tween  the  allowable  array  length,  which  causes  near-field  effects,  and 
array  flexure. 

NEAR-FIELD  OPERATION  AND  ARRAY  FLEXURE 

The  factors  that  produce  defocusing  are  (1)  near-field  operation, 
in  which  rays  to  or  from  the  earth  terminal  to  the  various  elements  of 
the  array  can  no  longer  be  considered  parallel,  and,  therefore,  of 
uniformly  progressive  length,  and  (2)  array  flexure,  which  introduces 
additional  path-length  changes.  Their  effect  is  analyzed  using  the 
geometry  illustrated  in  Fig.  3.3. 

The  earth,  of  radius  Rg,  is  centered  on  the  origin  0  of  the  x-y-z 
coordinate  system.  A  terminal  is  assumed  to  be  located  at  a  point  P, 
which  is  at  a  polar  angle  6  with  respect  to  the  z  axis  and  a  merid¬ 
ional  angle  <t>  with  respect  to  the  x-z  plane. 

The  PACSAT  array  is  assumed  to  consist  of  an  odd  number  N  of 
elements  with  the  center  element  at  the  origin  0'  of  the  x' -y* -z  co¬ 
ordinate  system.  The  center  element  is  at  an  altitude  h  above  the 
subsatellite  point  S.  The  x"  and  y'  axes  are  taken  parallel  to  the  x 
and  y  axes,  respectively,  and  the  z  and  z  axes  are  coincident.  The 
center  element  is  at  a  distance  from  the  earth  point  P,  so  the 
plane  containing  Iq  and  the  z  axis  is  also  at  a  meridional  angle  $ 
with  respect  to  the  x'-z'  plane.  This  path  from  P  to  O'  forms  a  nadir 
angle  e  with  respect  to  the  z'  axis. 

The  array  is  assumed  to  be  flexed  symmetrically  into  a  parabolic 
shape  and  its  axis  is  taken  as  the  tangent  at  its  center.  The  array 
axis  is  tilted  at  a  polar  angle  from  the  z  axis  and  is  at  a 
meridional  angle  v  with  respect  to  the  y'-z'  plane.  The  plane  con¬ 
taining  the  flexed  array  is  at  an  angle  id  with  respect  to  the  plane 
containing  the  array  and  z'  axes. 


I 

I 
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x 


The  kth  element  of  the  array  Is  at  the  point  K  a  distance  from 
the  center  element  and  a  distance  £rom  the  point  P.  The  lines 

from  0'  to  K  and  the  array  axis  form  an  angle  y  .  A  sphere  of  radius 
dK  is  centered  on  the  x'-y'-z'  coordinate  system;  its  intersections 
with  the  +z'  axis,  the  array  axis,  and  the  line  from  P  to  O'  are 
designated  by  T,  U,  and  V,  respectively. 

The  polar  angle  to  the  line  from  0'  to  K  together  with  <p 
and  y,  ,  which  have  already  been  defined,  form  the  sides  of  the  oblique 
spherical  triangle,  TUK;  the  included  angle  at  T  is  denoted  by 
The  great  circle  from  K  to  V  completes  another  oblique  spherical 
triangle,  TKV,  whose  sides  are  p^,  6^,  and  u-  e  and  whose  included 
angle  at  T  is  equal  toir/2  -  <j>  -  v  -  A^. 

The  geometry  shown  in  Fig.  3.3  is  general  and  bears  no  relation¬ 
ship  to  a  particular  orbit  for  PACSAT.  In  the  next  section,  when 
results  will  be  specialized  for  a  circular  equatorial  orbit,  the 
y-axis  will  become  the  earth's  polar  axis  and  the  x-z  plane  will 
become  the  equatorial  plane.  The  unusual  orientation  used  in  Fig.  3.3 
was  chosen  to  afford  the  clearest  perspective  of  the  various  angles 
and  distances  that  are  involved. 

In  the  following  analysis,  the  path  length  from  a  point  P^  to 

another  point  P£  via  the  kth  scatterer  will  be  calculated  as  a  power 

series  in  k,  as  was  indicated  in  Eq.  (3.2).  The  coefficient  C2  of  the 
o 

term  in  k  will  then  serve  to  express  the  defocusing  parameter  U  from 
Eq.  (3.15)  in  terms  of  the  significant  geometrical  factors  shown  in 
Fig.  3.3.  As  suggested  above,  the  subscripts  1  and  2  will  be  used  to 
identify  those  factors  that  are  associated  with  the  points  P^  and  P2, 
respectively.  The  expansion  will  include  near-field  effects  by  assum¬ 
ing  that  is  sufficiently  small  in  comparison  with  Iq  that  only  the 
leading  terms  in  d^Ag  need  be  retained.  Similarly,  flexure  effects 
will  be  included  by  assuming  that  is  sufficiently  small  that  only 
leading  terms  in  y^  need  be  retained. 

Basic  Expressions 

There  are  no  limitations  on  the  meridional  angle  <{>  of  the  point 
P.  Its  polar  angle  0,  however,  is  limited  to  the  value  at  which  P  is 
on  the  horizon.  This  occurs  when 
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cos  emaK  "  R  +  h 


(3*17) 


,  ..t,  o  tn  the  center  element  of 

the  triangle  00' P,  the  path  length  lQ 


the  array  la  given  by 


*2  „  r2  +  (R  +  h)2  -  2Re(Re  +  cos  6 
0  e  e 


(3.18) 


and  the  nadir  angle  e  by 


sin  e 


—  sin  8 


(3.19) 


The  path  length  from  P  to 
found  from  the  triangle  O'KP  as 


the  kth  element  of  the  array  Is 


«?  -  *5  *  4  -  “8 


(3.20) 


Expanding  In  powers  of  dk  leads  to 


c“*k  +  **0 


+  4  %  sin2  «k  +  0  4 


(3.21) 


The  ...  CO.  Sk  i.  «ochd  fee  the  .pherld  trt.ngle  VTh  .. 


cos  6k  -  -  cos  e  cos  Pk 


+  sin  e  sin  Pk  x 


+  w)  cos  A.  +  cos  (♦  +v>  sinXk]  (3*22> 


jsin  ($  +  v)  cos  +  cos  (<J> 
From  the  .pherlc.l  tri.ngle  T0K  it  1. 
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cos  *  cos  iji  cos  Yfc  -  sin  4)  sin  Y k  cos  w 
sin  y.  sin  w 

8ln  \  -  - sinTT  <3*23> 

k 

cos  Y^  *  cos  i}>  cos  +  sin  sin  and 

Substituting  Eqs.  (3,23)  in  Eq.  (3.22)  leads,  after  some  manipulation, 
to 


where 


cos  <$k  *  h  cos  Y^  +  sin  Yk  (3.24) 


h  »  -  cos  e  cos  <1/  +  sin  e  sin  (<p  +  v)  sin  ^ 


j  »  cos  e  sin  cos  u>  +  sin  e  sin  (<p  +  v)  cos  ^  cos  u>  (3.25) 


+  sin  t  cos  (<p  +  v)  sin  w 


Expanding  Eq.  (3.24)  in  powers  of  Yk  then  leads  to 


cos  -  h  +  j  Yk  -  j  h  y£  +  0  |y^| 


sin2  «k  -  (1  -  h2)  -  2hj  Yk  +  (h2  -  j2)  y2  +  0  (y*) 


(3.26) 


which  are  needed  in  Eq.  (3.21). 


Parameters  of  the  Flexed  Array 

The  array  of  length  L  Is  assumed  to  be  flexed  in  the  shape  of  a 
plane  parabola  of  maximum  deflection  D  as  shown  In  Fig.  3.4.  The  x 
axis  corresponds  to  the  array  axis  In  Fig.  3.3.  As  before,  the  number 
of  elements  N  Is  taken  as  odd  so  that  there  is  an  element  at  the  cen¬ 
ter.  The  array  length  Is  related  to  the  number  of  elements  and  their 
spacing  s  by 


L  «  (N  -  1 ) s 


(2.39) 


Array 

axis 


Fig.  3.4  —  Geometry  of  the  flexed  array 


To  facilitate  the  analysis,  it  will  be  assumed  that  the  array 
flexes  smoothly.  The  rectilinear  distance  sk  along  the  array  from  its 
center  to  the  kth  element,  which  is  given  by 


sk  -  sk,  k  -  -  2-=-^  ,  (3.27) 


can  then  be  written  as  the  curvilinear  distance 


s(t)  -  t, 


“  r 


(3.28) 


where  t  is  a  running  variable  along  the  array. 
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Let  L'  denote  the  length  of  the  array  projected  on  the  x  axis. 
It  can  be  shown  for  D  «  L  that  a  negligible  error  will  result  by 
using  L  for  V  .  The  array  is  then  described  by 


y(t)  -  x2(t) 

L 


(3.29) 


To  obtain  a  parametric  representation  of  a  point  on  the  array, 
note  that  the  differential  length  ds  along  the  array  can  be  written 


to  2 .  to 2  +  to 

\dt )  \dt  /  \dt  / 


(3.30) 


Using  Eq.  (3.28)  for  s  and  Eq.  (3.29)  for  y  yields 


(3.31) 


It  can  be  shown  that  the  solution  to  Eq.  (3.31)  is  given  by 


x(t)  -  t 


(3.32) 


Using  Eq.  (2.39)  to  eliminate  L  and  Eqs.  (3.27)  and  (3.28)  to  reintro¬ 
duce  k  leads  to 


„  _  Dt  Dk  „  . 

'7  «<»  -  i)  <(  1 


<3.33) 


Thus,  Eq.  (3.32)  can  be  written 
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\ 


sk  1 


|i32  +  0(^4) 


(3.34) 


and,  from  Eq.  (3.29), 


4  sk  Jb  1 1 


- 

3  ^ 


0(34) 


(3.35) 


The  distance  d^  and  the  angle  t0  ^e  kth  element  are  then  found  to 
be 


dk =  V  *k +  yk =  sk  [i  -  f +  °^4>]  (3*36> 

and 

7k  *  tan_1  *  4-4  |l  -  16  J>2  +  0(^4)|  (3.37) 

which,  noting  Eq.  (3.33),  are  the  desired  expansions  In  powers  of  k. 
The  Path-Length  Expansion 

2 

The  terms  cos  and  sin  6^  in  Eqs.  (3.26)  are  found  by  substi¬ 
tuting  from  Eq.  (3.37) 


cos  6k  =  h  +  4 -  8h  +  0C53) 
sin2  <$k  «  (1  -  h2)  -  8hjj4+  0(^2)  (3.38) 


These,  together  with  Eq.  (3.36)  for  d^,  then  give  the  path  length 
&k  In  Eq.  (3.21)  as 
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i,  =  £n  -  shk  + 
k  0 


2  2 
s  1  -  h 


4D 


l  2 

0  (N  -  1) 


hj.k3  +  0(k4) 


(N  -  l)2  lQ 


(3.39) 


where  Eq.  (3.33)  has  been  used  to  eliminate  3. 

3 

For  practical  values  of  D,  s,  N,  and  the  term  in  k  is  neg¬ 
ligible  so  the  total  path  length  Lk  from  the  point  P^  to  the  point 
via  the  kth  element  is  given  by 


L 


k 


-  s(h1  +  h2)k 


+ 


_ 4D 

(N  - 


1) 


2^1  +  ^2^ 


(3.40) 


where  the  subscripts  identify  terms  as  they  are  related  to  the  two 
points,  and  where  h  and  j  are  given  by  Eqs.  (3.25). 

A  comparison  of  Eq.  (3.40)  with  Eqs.  (3.2)  to  (3.5)  shows  that 


+ 


Cj  ■  »  -s(h1  +  h2) 


(3.41) 


c2 


4D 

(N  -  l)2 


<J1  +  j2) 


The  first  of  these  is  as  expected.  The  second  establishes  the  rela¬ 
tionships  between  the  basic  PACSAT  incidence  and  first-order  grating 
lobe  angles  a  and  0^,  respectively.  From  Eq.  (2.24)  it  is  seen  that 
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cos  a  =  -h^,  cos  6^  =  -t>2 


(3.42) 


The  Defocuslng  Parameter 

The  third  of  Eqs.  (3.41)  provides  the  coefficient  C2  needed  to 
calculate  the  defocuslng  parameter  U  from  Eq.  (3.15).  Although  this 
formulation,  which  is  expressed  in  terms  of  the  basic  geometric  system 
parameters,  is  suitable  for  numerical  evaluation,  it  does  not  lend 
itself  to  an  interpretation  of  how  the  near-field  and  flexure  param¬ 
eters  contribute  to  defocuslng.  To  obtain  such  a  form,  it  is  neces¬ 
sary  to  simplify  the  parameters  and  in  a  fashion  similar  to  that 
done  with  h^  and  h2  in  Eqs.  (3.42). 

The  simplified  geometry  of  the  tilted  and  flexed  array  is  shown 
in  Fig.  3.5,  which  is  similar  to  the  upper  portion  of  Fig.  3.3.  The 


Fig.  3.5  —  Simplified  geometry  of  the  tilted  and  flexed  array 
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flexed  array  is  shown  as  an  arc  from  the  origin  O'  passing  through  a 
point  K.  The  array  axis  and  the  flexed  array  then  define  the  plane  in 
which  the  array  is  tilted.  Note  that  the  angle  VO' U  is  equal  to 
n  -  a  or  n  -  fi ^  depending  on  whether  the  point  on  the  earth,  which  is 
not  shown,  is  the  transmitter  or  the  receiver,  respectively.  Thus, 
the  array  axis  and  t^e  line  of  length  Iq  from  O'  to  the  earth  point 
define  the  plane  of  illumination  (or  scattering). 

Let  the  c;  denote  the  angle  between  the  plane  of  illumination  (or 
scattering)  and  the  plane  of  flexure.  It  can  be  shown  that 

»  sin  a  cos  ^2  =  sin  0^  cos  4 ^  (3.43) 

Then,  using  Eq .  (2.39)  to  reintroduce  the  array  length  L  and  Eqs. 
(3.42)  and  (3.43)  to  eliminate  h  and  j  brings  the  defocusing  parameter 
in  Eq.  (3.15)  into  the  form 


U2  - 


.  2 
sin  a 


Sin2  0, 

~2STi 


^3.  (sin  a  cos  +  sin  0^  cos  C 


(3.44) 


where  the  third  of  Eqs.  (3.41)  was  used  for  C2 • 

The  first  group  of  terms  in  Eq.  (3.44)  exhibits  the  effects  of 

near-field  operation.  As  the  array  length  L  is  increased  in  relation 

to  the  wavelength  A  and  the  array  center  path  lengths  and  , 

I  t>2 

these  terms  increase  the  value  of  the  defoc  sing  parameter.  It  is 

interesting  to  note  that  L  sin  a  and  L  sin  0 ^  are  the  projected 

lengths  of  the  array  in  the  directions  of  illumination  and  scattering, 

respectively.  Thus,  near-field  effects  are  reduced  as  end-fire 

operation  is  approached.  One  of  the  adverse  effects  of  libration  is 

that  the  near-field  effect  is  exacerbated  when  the  array  is  at  the 

less  favorable  orientations.  It  should  be  noted  that  a  and  0^  vary 

during  the  transit  of  an  array  in  orbit,  even  if  there  is  no  libration. 


i 


}  2 

For  equal  path  lengths,  the  term  {(L  sin  a)4  +  (L  sin  8^)  ]/2A^ 
is  equal  to  the  average  of  the  ratios  of  the  squares  of  the  projected 
array  length  to  the  operating  wavelength.  Thus,  the  near-field  term 
in  Eq.  (3.44),  when  equated  to  unity,  is  analogous  to  the  approxima¬ 
tion  R  =  D 2/X  used  in  aperture  theory  to  denote  the  range  R  at  which 
the  transition  from  the  Fresnel  (near-field)  to  the  Fraunhofer  (far- 
field)  regions  occurs.  That  is  to  say,  the  value  U  =  1.0  mentioned 
above  as  the  desirable  maximum  allowable  value  for  U  is  consistent,  in 
the  absence  of  flexure,  to  operation  of  the  PACSAT  array  in  the  tran¬ 
sition  region. 

The  second  group  of  terms  in  Eq.  (3.44)  displays  the  effects  of 
array  flexure.  Unlike  the  near-field  terms,  which  are  nonnegative, 
the  flexure  terms  can  be  negative  and  can,  therefore,  actually  compen¬ 
sate  for  near-field  effects  to  some  extent.  Thus,  they  can  focus  the 
array  as  well  as  defocus  it.  To  appreciate  this,  note  from  Fig.  3.5 
that  if  4  is  near  zero,  the  array  is  flexed  so  as  to  present  a  con¬ 
cave,  or  focusing,  aspect  along  the  line  of  sight.  Conversely,  if 
4  is  near  x,  the  aspect  is  convex,  or  defocusing.  As  with  the  near¬ 
field  terms,  the  factors  sin  a  and  sin  8^  indicate  that  flexure  ef¬ 
fects  are  reduced  as  end-fire  operation  is  approached. 

The  Flexure-Length  Tradeoff 

The  form  of  the  defocusing  parameter  in  Eq.  (3.44)  suggests  that 
a  tradeoff  exists  between  the  defocusing  produced  by  flexure  and  that 
produced  by  near-field  operation.  To  illustrate  this  tradeoff,  con¬ 
sider  a  pair  of  users  at  38  deg  latitude  separated  by  38  deg  in  longi¬ 
tude.  (These  locations  are  3310  km  apart  at  a  mid-U.S.  latitude.  They 
are  roughly  midway  among  the  various  user  pairs  in  the  next  section 
(see  Fig.  4.5)  that  illustrate  various  system  parameters  for  PACSAT  in 
a  circular  equatorial  orbit.) 

Let  PACSAT  be  at  an  altitude  h  above  the  equator  at  a  longitude 
midway  between  the  users.  Further,  let  the  array  be  tilted,  i.e., 
librated,  toward  the  north  at  angle  iJj  so  that  it  presents  a  maximum 
projected  length  toward  the  users  (v  *  0  in  Fig.  3.5).  Finally,  let 
the  array  be  flexed  toward  the  south  so  as  to  present  a  convex  or 
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defocusing  aspect  toward  the  users  (u  =  r  In  Fig.  3.5).  Although  this 
is  probably  not  a  worst-case  configuration  (the  users  should  be  nearer 
one  another  and  at  some  pessimum  latitude),  it  is  nonetheless  one  that 
is  highly  unfavorable  and,  therefore,  suitable  for  illustrating  the 
adverse  effects  of  defocusing. 

The  tradeoff  is  plotted  in  Fig.  3.6,  which  shows  how  the  nor- 
2 

malized  array  length  L  /\^Re  can  be  traded  off  against  the  normalized 
flexure  D/A^  for  a  constant  value  of  1.0  for  the  defocusing  parameter. 

The  most  striking  thing  to  be  noted  from  Fig.  3.6  is  that  the 
flexure  cannot  be  allowed  to  exceed  a  fraction  of  a  wavelength  for  the 
range  of  parameters  considered,  even  if  there  were  no  near-field  ef¬ 
fects.  If  practical  array  lengths  are  assumed,  the  near-field  effect 
further  reduces  the  flexure  that  can  be  tolerated.  The  straightness 
criterion  is  therefore  seen  to  be  stringent;  for  example,  an  eighth  of 
a  wavelength  is  less  than  5  mm  at  8  GHz.  Orbiting  an  array  that  de¬ 
parts  from  linearity  by  no  more  than  such  an  amount  over  a  length  of 
1000  m  or  more  will  be  a  formidable  task. 


Fig.  3.6  —  Near-worst-case,  flexure/length  tradeoff:  users  at  38  deg  latitude 
separated  38  deg  in  longitude  and  array  at  midlongitude  over  equator 


Normalized  array  length  L2/AR 
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The  significance  of  the  limitation  in  the  normalized  array  length 

2 

L  A^R  can  be  appreciated  from  Fig.  3.7,  which  shows  how  it  is 


Approximate  scattering  cross  section,  m2 
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Array  length  L,  m 

Fig,  3.7  —  Normalized  array  length 


-50- 


re  la  ted  to  the  actual  array  length  at  a  given  operating  frequency. 
Consider,  for  instance,  a  system  operating  at  an  orbital  altitude  of 
8000  km.  If  it  were  felt  that  the  array  flexure  could  be  limited  to 
an  eighth  of  a  wavelength  and  the  array  libration  to  two  degrees,  the 
tradeoff  in  Fig.  3.6  shows  that  the  normalized  array  length  must  not 
exceed  a  value  of  about  5.2.  From  Fig.  3.7  it  is  seen  that  the  array 
length  would  be  limited  to  about  1100  m  if  the  operating  frequency 

were  8  GHz.  The  scale  at  the  top  of  Fig.  3.7,  which  is  constructed 

2 

according  to  the  approximation  a  =  0.6L  (see  Eq.  (4.15))  shows  that 
an  array  of  that  length  would  have  a  scattering  cross  section  of  about 
725,000  m2. 

The  specific  values  cited  above  must,  of  course,  be  taken  in 
context  because  they  relate  to  a  specific  user  and  orbital  geometry 
that  may  or  may  not  be  worst  case;  to  an  array  configuration  and  be¬ 
havior  (i.e.,  X/8  flexure  and  2  deg  libration)  that  may  or  may  not  be 
achieved;  and  to  a  defocusing  criterion  (i.e.,  U  =  1.0)  that  is  sub¬ 
jective  to  some  extent.  The  important  conclusions  are  that  (1)  the 
array  must  be  straight  to  within  a  fraction  of  a  wavelength;  (2)  the 

libration  should  be  limited  to  a  few  degrees,  at  most;*  and  (3)  the 

2 

normalized  array  length  L  /^Re  must  not  exceed  a  value  on  the  order 
of  10.  If  these  conditions  are  not  satisfied,  significant  reductions 
in  scattering  cross  section  and  distortions  of  the  scattered  beam  (in 
the  vicinity  of  the  mainlobe)  will  occur  under  some  unfavorable  cir¬ 
cumstances  . 

*It  will  be  shown  in  the  next  section  that  excessive  libration 
also  exacerbates  the  system  bandwidth  requirement  problem. 
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IV.  PACSAT  IN  A  CIRCULAR  EQUATORIAL  ORBIT 


The  desirability  of  a  circular  equatorial  orbit  for  PACSAT  was 
mentioned  in  the  Synopsis.  In  fact,  circularity  is  mandatory  because 
gravity-gradient  stabilization,  which  is  used  to  keep  the  array  erect, 
will  not  work  in  an  elliptical  orbit.  Placing  the  circular  orbit  in 
the  equatorial  plane  offers  important  system  advantages  that  outweigh 
the  loss  of  polar  coverage  provided  only  by  satellites  in  inclined 
orbits;  at  the  orbital  altitudes  suitable  for  use  by  PACSAT,  coverage 
can  be  provided  up  to  the  midlatitudes.  Thus,  strategic  command  and 
control  communications  that  require  CONUS  coverage  can  be  satisfied. 

A  practical  system  must  incorporate  enough  satellites  to  provide 
continuous  service.  To  form  a  communication  link,  a  pair  of  users 
must  acquire  a  predetermined  satellite  as  it  comes  into  mutual  view, 
then  track  it  during  its  transit.  At  some  point  the  satellite  must  be 
relinquished  and  another  acquired  to  maintain  continuity.  Inasmuch  as 
satellites  in  circular  equatorial  orbits  make  identical  transits  as 
viewed  by  a  given  user,  the  antenna  pointing  problem  for  acquisition 
and  tracking  is  simplified.  Furthermore,  only  a  limited  segment  of 
the  celestial  sphere  need  be  viewed,  th>.  reby  facilitating  terminal 
siting  and  antenna  design  in  some  instances. 

Another  advantage  of  the  circu  ar  equatorial  orbit  is  that  the 
PACSAT  design  uses  a  frequency  steerable  beam  to  form  a  communication 
link.  As  noted  in  connection  with  Eq.  (1.1),  the  ap. vopriate  operat¬ 
ing  wavelength  depends  on  the  illumination  and  scattering  angles  a 
and  3^.  If  these  angles  are  allowed  to  assume  the  wide  range  of 
values  that  can  result  when  inclined  orbits  are  used,  the  operating 
bandwidth  can  greatly  exceed  available  allocations  It  is  a  charac¬ 
teristic  of  the  circular  equatorial  orbit  that,  for  users  at  midlati¬ 
tudes,  the  values  of  a  and  3^  vary  in  such  a  way  during  the  course  of 
a  satellite  transit  that  the  sum  of  their  cosines  is  nearly  constant. 
(In  fact,  the  bandwidth  requirement  will  depend  principally  on  the 
amount  of  libration  the  array  experiences.)  Thus,  the  use  of  circular 
equatorial  orbits  reduces  the  required  r.parating  bandwidth  to  a 
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practical  amount  and  eases  the  frequency,  as  well  as  the  spatial, 
acquisition  and  tracking  problem. 

In  this  section,  the  orbital  properties  of  satellites  in  circular 
equatorial  orbits  are  reviewed  and  the  important  geometrical  relation¬ 
ships  to  viewers  at  typical  locations  are  developed.  This  is  followed 
by  an  examination  of  the  pertinent  electromagnetic  or  communication 
characteristics.  These  are  then  used  to  illustrate  a  typical  array 
design.  Finally,  the  link  equation  is  derived  and  typical  power  bud¬ 
gets  are  presented. 

GEOMETRICAL  AND  ORBITAL  PROPERTIES 


Orbital  Rate  and  Period 

The  links  for  PACSAT  in  a  circular  equatorial  orbit  are  illus¬ 
trated  in  Fig.  4.1.  The  earth  is  centered  on  the  origin  of  the  x-y-z 
coordinate  system  with  the  polar  axis  coinciding  with  the  z  axis.  The 
earth  is  rotating  at  an  angular  rate  SI  equal  to  7.2921  x  10-^  rad/sec, 
which  is  the  rate  corresponding  to  the  23  hr  56  min  4.1  sec  sidereal 
day.  However,  the  coordinate  system  is  assumed  to  be  fixed  to  the 
earth  with  the  Greenwich  meridian  in  the  x-z  plane.  PACSAT  shown 
in  the  x-y,  or  equatorial,  plane  at  a  constant  altitude  h  the 

surface  of  the  earth  and  is  assumed  to  be  in  a  progressive  orbit, 
i.e.,  rotating  in  the  same  direction  as  the  earth. 

The  orbital  angular  rate  for  an  earth  satellite  relative  to  iner¬ 
tial  space  is  given  by 


Co 


rad/ sec 


(4.1) 


where  r  is  the  orbital  radius  and  p  =  1.4076  x  10^  ft^/sec^.  Con¬ 
verting  to  kilometers  and  introducing  the  altitude  h  and  the  radius  of 
the  earth  Rg  yields 


6.3135  x  10 


(h  +  R  ) 
e 


3/2 


U) 


rad/sec 


(4.2) 


PACSAT 


Fig.  4.1  —  Link  geometry  for  PACSAT  in  a  circular  equatorial  orbit 


where  Rg  *  6378  km.  This  orbital  period  2tt /u>  in  hours  is  plotted  on 
Fig.  4.2  for  a  range  of  orbital  altitudes  from  3000  to  10,000  km, 
which  encompasses  those  of  interest  for  PACSAT.  Let  Qg  denote  the 
longitude  of  the  satellite.  Then,  the  orbital  rate  relative  to  the 


3000  4000  5000  6000  7000  8000  9000  10,000 

Orbital  altitude  h,  km 


Fig.  4.2  —  Orbital  period  for  a  satellite  in  a  circular  orbit 


earth  Is  given  by 

<3S  =  u  ~  0  (4.3) 

This  is  plotted  on  Fig.  4.3  in  deg/min. 

Coverage 

Typical  users  are  shown  in  Fig.  4.1  at  points  and  ?2,  which 
are  at  latitudes  and  longitudes  0^  and  * ©£*  respectively.  The 
lines  of  sight  from  these  users  to  the  array  form  a  bistatic  scattering 
angle  y  and  the  bearing  angles  e  and  z  .  If  the  array  is  erect, 
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Orbital  altitude  h,  km 


Fig.  4.3  —  Orbital  angular  rate  relative  to  the  earth  for  a 
satellite  in  a  circular  equatorial  orbit 


these  become  the  incidence  and  scattering  angles  a  and  3^  as  shown  in 
Fig.  1.1.  The  great  circles  from  these  points  to  the  subsatellite 
point  form  polar  and  meridional  angles  6^,  and  6^,  respec¬ 
tively,  as  also  shown  in  Fig.  3.3.  These  angles  are  related  to  the 
various  latitudes  and  longitudes  by 


COS  0  “  COS  (0p  -  0g)  COS  *p 


tan 


tan  $p 

sin  (0 p  -  Og) 


(4.4) 


where  the  appropriate  numerical  subscripts  are  to  be  added  as  required. 
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The  satellite  appears  at  the  horizon  to  a  given  user  when  its 
line  of  sight  is  tangent  to  the  surface  of  the  earth.  At  that  point 


cos  9 


max 


R 

e 

R  +  h 
e 


(3.17) 


so,  using  Eq.  (4.4),  the  satellite  longitudes  at  the  extremes  of  its 
visibility  are  given  by 


0g  =  0p  ±  cos 


R 

_ e _ 

(Rg  +  h)  cos  $p) 


(4.5) 


For  the  progressive  satellite  motion  depicted  in  Fig.  4.1,  the  longi¬ 
tude  at  which  the  satellite  rises  is  obtained  by  using  the  minus  sign 
in  Eq.  (4.5);  the  longitude  of  setting  then  comes  from  using  the  plus 
sign. 

The  coverage  afforded  by  a  satellite  in  an  equatorial  orbit  is 
shown  in  Fig.  4.4,  which  is  a  portion  of  a  polar  stereographic  projec¬ 
tion  [6].  The  plane  of  the  projection  is  tangent  to  the  earth  at  one 
pole,  and  the  point  from  which  the  projection  is  developed  is  the 
other  pole.  Such  a  projection  is  conformal  and,  hence,  has  the  prop¬ 
erty  that  circles  on  the  earth  are  projected  into  circles.  Thus, 
coverage  areas,  which,  of  course,  are  circles,  are  projected  as 
circles  rather  than  the  distorted  figures  usually  seen  on  other 
projections. 

The  projection  in  Fig.  4.4  is  limited  to  one  hemisphere,  so  the 
outer  circular  arc  of  the  plot  is  the  equator.  The  coverage  circles 
show  the  portion  of  the  earth  visible  at  various  altitudes  from  a 
satellite  at  the  indicated  subsatellite  longitude.  With  respect  to 
the  coverage  circle  associated  with  a  particular  orbital  altitude,  a 
user  will  be  in  view  of  the  satellite  if  he  is  within  the  coverage 
circle.  If  the  user  is  on  the  coverage  circle,  the  satellite  will 
appear  to  be  on  the  horizon.  The  two  longitudes  at  which  a  user  at  a 
given  latitude  can  be  on  the  coverage  circle  are  separated  by  the 
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•  Points  used  in  the  text  example. 

Fig.  4.4  —  Coverage  contours  for  a  satellite  in  an  equatorial  orbit 


difference  In  rising  and  setting  satellite  longitudes  given  by  Eq. 
(4.5).  For  example,  consider  a  user  at  a  latitude  of  40  deg  and  a 
satellite  at  an  altitude  of  6000  km;  the  two  relative  satellite  longi¬ 
tudes  at  which  that  user  would  observe  the  satellite  at  the  horizon 
are  shown  by  dots  on  Fig.  4.4.  The  separation  in  longitude  is  seen  to 
be  95.5  deg,  which  is  equal  to  twice  the  inverse  cosine  term  in  Eq. 
(4.5). 

Because  of  the  potential  application  of  PACSAT  within  the  con¬ 
tinental  United  States,  it  is  convenient  to  define  a  generalized  ser¬ 
vice  area  that  approximates  the  United  States  in  shape.  User  loca¬ 
tions  at  the  extremities  of  this  area  can  then  be  paired  to  illustrate 
the  significant  system  parameters  to  be  developed  below.  The  service 
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area  will  be  assumed  to  extend  45  deg  in  longitude  between  30  and  45 
deg  N  latitude.  This  area  is  shown  superposed  over  the  United  States 
on  the  polar  stereographic  projection  in  Fig.  4.5.  The  circles  at  the 
corners  of  this  area  are  user  locations  that  form  the  following  pairs: 

1.  Both  users  at  30  deg  latitude  separated  by  45  deg  in  longi¬ 
tude. 

2.  One  user  at  30  deg  latitude  and  the  other  at  45  deg  latitude 
separated  by  45  deg  in  longitude. 

3.  Both  users  at  45  deg  latitude  separated  by  45  deg  in  longi¬ 
tude  . 

These  will  be  used  to  illustrate  some  of  the  system  parameters  de¬ 
veloped  in  the  following.  The  results  that  do  not  require  specifi¬ 
cation  of  a  user  pair  will  be  presented  for  users  at  30  deg  and  45  deg 
latitude.  The  user  pair  located  at  38  deg  latitude  separated  by  38 
deg  in  longitude  chosen  to  illustrate  the  flexure/length  tradeoff  in 
Fig.  3.6  is  indicated  by  stars  in  Fig.  4.5. 

Transit  Trajectories 

The  line  of  sight  from  a  point  on  the  earth  to  a  satellite  can  be 
characterized  by  its  azimuth  A  and  elevation  E.  These  are  defined  in 
Fig.  4.1. 

The  azimuth  is  seen  to  be  independent  of  the  satellite  altitude 
and  is  given  by 


cot  (it  -  A)  *  sin  Gp  cot  (Qp  -  0g)  (4.6) 

in  the  northern  hemisphere  and  by 

cot  A  *  -  sin  cot  (®p_0s)  (4.7) 


in  the  southern  hemisphere.  Although  these  forms  are  essentially 
identical,  they  are  listed  separately  to  Insure  that  using  principal 
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values  results  In  azimuths  that  are  consistent  with  the  definition. 

The  elevation  is  given  by 

E  =  j  -  9  -  e  (4.8) 

where  6  is  given  by  Eq.  (4.4)  and  e  by  Eq.  (3.19). 

Azimuth  and  elevation  trajectories  are  plotted  in  Figs.  4.6a  and 
4.6b  for  users  at  latitudes  of  30  and  45  deg  north  for  the  range  of 
orbital  altitudes  from  3000  to  10,000  km.  In  general,  for  users  in  the 
northern  hemisphere,  the  satellite  rises  from  the  south  of  west,  reaches 
its  maximum  elevation  when  due  south,  then  sets  to  the  south  of  east. 


Relative  satellite  longitude  0,-©|,  deg 


Fig  4.6a  —  Azimuth  and  elevation  trajectories  for  a  satellite  in  a  circular 
equatorial  orbit  as  viewed  by  a  user  at  30  deg  N  latitude 
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Fig.  4.6b  —  Azimuth  and  elevation  trajectories  for  a  satellite  in  a  circular 
equatorial  orbit  as  viewed  by  a  user  at  45  deg  N  latitude 

Mutual  Visibility 

The  coverage  of  Interest  for  communication  purposes  is  that 
within  which  the  satellite  is  in  view  by  a  pair  of  users.  Figure  4.4 
can  be  used  to  determine  if  the  users  have  mutual  view  of  a  satellite 
by  seeing  if  they  are  both  within  the  coverage  area.  Because  of  the 
relative  motion  of  the  satellite,  the  coverage  area  can  be  thought  to 
rotate  in  a  counterclockwise  direction,  i.e.,  toward  the  east,  as  the 
points  denoting  the  users  remain  fixed  on  the  projection  grid.  Con¬ 
versely,  the  coverage  area  can  be  thought  to  remain  fixed  as  the  users 
rotate  in  a  clockwise,  or  westerly,  direction.  By  either  method,  the 


Azimuth  A,  deg  (referenced  to  north) 
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relative  longitude  at  which  the  satellite  first  appears  in  view  of 
both  users  and  then  passes  from  view  of  both  can  be  estimated. 

The  pairs  of  rising  and  setting  relative  satellite  longitudes  can 
also  be  calculated  from  Eq.  (4.5)  and  the  pair  representing  the  ex¬ 
tremes  of  mutual  visibility  selected  by  inspection.  It  should  be 
noted  that  the  satellite  does  not  necessarily  appear  to  either  rise  or 
set  first  as  seen  by  the  more  westerly  user.  The  more  easterly  user 
may  observe  the  satellite  to  rise  first  if  he  is  sufficiently  nearer 
the  equator;  similarly,  he  may  observe  it  to  set  first  if  he  is  suf¬ 
ficiently  nearer  the  pole. 

The  arc  of  mutual  visibility,  i.e.,  the  span  of  satellite  longi¬ 
tudes  over  which  the  satellite  is  visible  to  both  users,  is  shown  in 
Figs.  4.7a  to  4.7c  for  the  user  pairs  designated  above.  As  might  be 


Orbital  altitude  h.  km 


Fig.  4.7a  -  Arc  of  mutual  visibility  for  a  satellite  in  a  circular 
equatorial  orbit  as  viewed  by  users  both  at  30  deg  latitude 


Orbital  arc  of  mutual  visibility,  deg  Orbital  arc  of  mutual  visibility,  deg 


expected,  the  arc  of  mutual  visibility  is  greater  for  the  users  nearer 
the  equator  and  increases  as  the  satellite  orbital  altitude  is  raised 
or  as  the  users  are  brought  closer  in  longitude.  Plots  such  as  these 
are  useful  in  estimating  the  minimum  number  of  satellites  required  to 
maintain  continuous  coverage. 

The  mutual  viewing  time  is  obtained  by  applying  the  apparent 
orbital  angular  rate  from  Eq.  (4.3)  to  the  arcs  of  mutual  visibil¬ 
ity.  These  are  presented  in  Figs.  4.8a  to  4.8c  for  the  same  user 
pairs  as  above. 

For  user  pairs  at  other  locations,  the  arc  of  mutual  visibility 
can  be  inferred  from  Fig.  4.4  by  the  procedure  previously  outlined. 

The  mutual  viewing  time  can  then  be  calculated  using  the  relative 
orbital  angular  rate  from  Fig.  4.3. 


Fig.  4.8a  —  Mutual  viewing  time  for  a  satellite  in  a  circular 
equatorial  orbit  as  viewed  by  users  both  at  30  deg  latitude 


Mutual  viewing  time,  mm  Mutual  viewing  time,  mm 
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Orbital  altitude  h,  km 


ig.  4.8b  -  Mutual  viewing  time  for  a  satellite  in  a  circular 
equatorial  orbit  as  viewed  by  one  user  at  30  rinn 
latitude  and  the  other  at  45  deg  latitude 


Fig.  4.8c  —  Mutual  viewing  time  for  a  satellite  in  a  circular 
equatorial  orbit  as  viewed  by  users  both  at  45  deg  latitude 
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ARRAY  OPERATING  PARAMETERS 


Bistatic  Scattering  Angle 

The  relevance  of  the  blstatic  scattering  angle  to  the  scattering 
cross  section  of  the  array  was  established  In  Sec.  II.  It  can  be 
shown,  from  Fig.  4.1,  to  be  given  by 


cos  Y 


+  Jt2  -  4  R2  sin2  (T/2) 
C>2  e 


2  i. 


0 „ 


(4.9) 


where  the  line-of-sight  path  lengths  and  to  the  center  of 

U1  u2 

the  array  are  given  by  Eq.  (3.18).  The  angle  T,  not  shown  in  Fig. 
4.1,  is  the  central  angle  between  the  users  and  is  given  by 


Satellite  longitude  relative  to  mid-user  longitude,  deg 


Fig.  4.9a  -  Bistatic  scattering  angle  for  a  satellite  in  a  circular 
equatorial  orbit  as  viewed  by  users  both  at  30  deg  latitude 
and  separated  by  45  deg  in  longitude 


Satellite  longitude  relative  to  mid-user  longitude,  deg 


Fig.  4.9c  -  Bistatic  scattering  angle  for  a  satellite  in  a  circular 
equatorial  orbit  as  viewed  by  users  both  at  45  deg  latitude 
and  separated  by  45  deg  in  longitude 
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cos  T  =  sin  4  ^  sin  $2  +  cos  cos  ^2  cos  ^2  ~  (4.10) 

The  bistatic  scattering  angle  is  shown  in  Figs.  4.9a  to  4.9c  as  a 
function  of  satellite  longitude  relative  to  the  midlongitude  of  a  pair 
of  users  for  a  range  of  satellite  orbital  altitudes.  The  plots  are 
for  the  three  cases  of  the  user  pairs  at  the  extremities  of  the  gen¬ 
eralized  service  area  depicted  in  Fig.  4.5.  When  the  users  are  at  the 
same  latitude,  the  maximum  value  of  the  bistatic  scattering  angle  is 
reached  when  the  satellite  is  midway  between  the  users. 


It  was  shown  in  Eq.  (3.41)  that  the  operating  wavelength  A^  re¬ 
quired  to  form  a  link  between  a  pair  of  users  is  given  by 


XL  =  -s(h^  +  h2) 


(3.41) 


where  s  is  the  element  spacing  and  the  h’s  are  given  by  the  first  of 
Eqs .  (3.25).  The  quantities  h^  and  h£  are  expressed  in  terms  of  the 
relevant  geometrical  factors  illustrated  in  Fig.  3.3.  The  latitudes 
and  longitudes  of  the  users  and  the  satellite  can  be  introduced  by 
using  Eqs.  (4.4),  if  desired.  The  end-fire  operating  wavelength  A^, 
in  Eq.  (2.39),  can  be  used  to  eliminate  s  yielding 


-  V2  (^  +  h2) 


(3.42) 


The  relationship  expressed  by  Eq.  (3.42)  is  purely  geometric  and 
does  not  include  the  dynamics  of  satellite  and  user  motion.  Thus,  the 
equation  is  correct  only  when  applied  to  the  case  of  a  geostationary, 
nonlibrating  satellite  as  viewed  by  fixed  users.  In  that  event,  the 
operating  wavelength  A^  given  by  Eq.  (3.42)  is  not  only  correct  but  is 
also  constant  and  the  same  in  both  directions  of  propagation. 
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The  Doppler  shift  introduced  by  user  motion  is  familiar  and 
easily  accommodated.  The  effect  of  orbital  motion  is  complex  and  will 
not  be  considered  here;  a  detailed  analysis  will  appear  in  a  forth¬ 
coming  report  where  its  effects  on  signal  acquisition  and  tracking  are 
considered  [7].*  Fortunately,  the  resulting  difference  in  operating 
wavelength  is  small  for  the  orbital  altitudes  of  interest  so  that  Eq. 
(3.42)  is  adequate  for  the  purpose  of  estimating  the  total  operating 
bandwidth  that  may  be  required. 

The  extremes  of  relative  operating  wavelength  X  ^ /X ^  and  its  equiv¬ 
alent  f^/fg  have  been  calculated  using  Eq.  (3.42)  and  are  shown  in 
Figs.  4.10a  to  4.10c  as  a  function  of  libration  angle  for  a  range  of 
satellite  orbital  altitudes.  The  computational  procedure  was  to 
calculate  the  maximum  and  minimum  values  of  that  were  experi¬ 

enced  during  a  fictitious  satellite  transit  in  which  the  tilt  angle  i|; 
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Fig.  4.10a  —  Extremes  of  operating  wavelength  and  frequency  for  a  librating 
PACSAT  in  a  circular  equatorial  orbit  as  viewed  by  users  both  at  30  deg 
latitude  and  separated  by  45  deg  in  longitude 


The  major  consequence  is  that  the  path  is  no  longer  bilateral. 
Also,  there  Is  a  Doppler  smear,  aggravated  by  libration,  that  develops 
because  of  differential  Doppler  shifts  along  the  array. 
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Fig.  4.10b  —  Extremes  of  operating  wavelength  and  frequency  for  a  libratir 
PACSAT  in  a  circular  equatorial  orbit  as  viewed  by  one  user  at  30  deg 
latitude  and  the  other  at  45  deg  latitude,  separated  by  45  deg  in  longitudt 
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Fig.  4.10c  —  Extremes  of  operating  wavelength  and  frequency  for  a  librating 
PACSAT  in  a  circular  equatorial  orbit  as  viewed  by  users  both  at  45  deg 
latitude  and  separated  by  45  deg  in  longitude 


aa  Ratio  of  operating  to  end-fire 

Ratio  of  operating  to  end-fire  frequencies  fi/fo 

frequencies  fi/fo 


was  fixed  in  a  given  plane  u).  This  was  repeated  over  the  360  deg 
range  of  w  to  determine  the  extrema.  (The  existence  of  a  spread  in 
operating  wavelength  even  when  there  is  no  tilt  merely  reflects  the 
changing  geometry  characteristic  of  the  basic  satellite  transit.) 

This  procedure  obviously  ignores  the  physics  of  libration.  It 
simply  determines  the  extremes  of  operating  wavelength  that  might  be 
experienced  during  the  transit  of  a  satellite  with  a  given  tilt  angle 
i if  it  were  allowed  to  assume  all  possible  orientations  m.  The  actual 
libration  of  a  rigid  body  is  much  more  complex.  The  motion  can  be  re¬ 
solved  into  in-plane  and  out-of-plane  oscillations  whose  periods  are 
equal,  respectively,  to  1/^3  and  one-half  the  orbital  period.  The  re¬ 
lative  amplitudes  and  phases  are  different,  in  general,  and  are  coupled 
in  a  complex  manner.  It  suffices  here  to  use  the  simple  worst-case 
procedure  to  determine  the  maximum  operating  bandwidth  required  for  a 
given  tilt  angle,  which  is  loosely  referred  to  as  the  libration  angle. 

The  phenomenon  of  nearly  constant  operating  wavelength  alluded  to 
in  the  introduction  to  this  section  is  apparent  from  the  no-libration 
cases  in  which  the  minima  and  maxima  experienced  during  the  course  of 
a  transit  usually  differ  by  much  less  than  one  percent.  Introducing 
libration  greatly  increases  the  range  of  wavelengths  that  might  be 
required.  If,  in  addition,  complete  freedom  of  user  locations  within 
the  satellite  field  of  view  were  allowed,  bandwidths  of  20  to  30  per¬ 
cent  or  more  would  be  required.  A  typical  frequency  allocation  suit- 
able  for  use  by  a  passive  satellite  system  is  only  6  percent  wide,  so 
it  is  apparent  that  libration  must  be  limited  to  only  a  few  degrees  at 
most  and  that  some  user  geometries  must  be  avoided.  (The  latter  may 
not  prove  to  be  a  serious  limitation  as  it  amounts  to  avoiding  close 
user  spacings  and  equatorial  latitudes.)  The  higher  orbital  altitudes 
are  clearly  preferable  from  this  point  of  view. 

Typical  Array  Design 

An  expression  for  the  scattering  cross  section  of  the  PACSAT 
array  in  the  far  field  was  developed  in  Sec.  II  and  is  given  by  Eq. 

*7900  MHz  to  84  MHz,  earth-to-space,  for  example. 
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(2.37).  Using  Eq.  (2.40),  it  can  be  brought  into  the  form 


o 

L2(a/A0) 


it  a 


(4.11) 


which  is  plotted  in  Fig.  2.6.  The  normalized  scattering  cross  section 
2 

of  a  sphere  o  /  it  a  is  obtained  from  Eq.  (2.18)  and  is  a  function  of 
LK 

the  bistatic  scattering  angle,  denoted  by  9  in  that  section  and  by  y 
in  this,  and  the  dimensionless  parameter 

V  -■  Hr  (4*12) 

where  A  is  the  operating  wavelength  (see  Eq.  (2.2))  and  a  is  the 

2 

sphere  radius.  For  computational  purposes,  aLR/na  is  given  by  the 

polynomial  Eq.  (A.l)  in  the  appendix. 

A  suitable  form  of  the  array  scattering  cross  section  for  design 
2 

purposes  is  o/L  ,  where  L  is  the  array  length.  To  arrive  at  this 
2 

from  o/L  (a/A^),  it  is  necessary  to  specify  the  sphere  radius  a  and 
the  end-fire  wavelength  A^.  The  latter  is  related  to  the  sphere  spac¬ 
ing  s  by 


s  -  AQ/2  (2.40) 

2 

Thus,  if  a  and  s  (and  hence  A^)  were  specified,  o/L  could  be  calcu¬ 
lated  as  a  function  of  A  and  0  by  bringing  Eq.  (4.11)  into  the  form 


o 


1  /v  %2 

2 7  <k0a) 


(4.13) 


A  typical  design  procedure  for  specifying  a  and  s  will  be  illus¬ 
trated  for  a  nonlibrating  PACSAT  intended  to  serve  the  area  shown  in 
Fig.  4.5.  It  is  based  on  the  desire  to  keep  the  system  operating 
point,  as  determined  by  A  and  0,  as  near  as  possible  to  the  ridge 
marked  on  the  surface  of  Fig.  2.6.  The  bistatic  scattering  angle  0 
is,  of  course,  fixed  by  the  satellite-user  geometry  so  that  only  the 
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normalized  sphere  size  kga  is  available  f : ..  specification.  Let  the 
subscript,  des,  denote  a  value  at  a  design  point.  Then 


kQa 


2na 

X 


2na  des 

X  .  Xo 
des  0 


'0 

X 


(kna)  (X/X  ) 

des  des 


T7x 


(4.14) 


0 


The  design  value  (^oa^des  he  selected  by  noting,  from  Fig. 

2.6,  that  the  array  scattering  cross  section  will  be  more  or  less 
optimized  for  lesser  values  jf  0  if  it  lies  on  the  ridge  when  0 
attains  its  maximum  value.  For  users  within  the  designated  service 
area,  this  maximum  occurs  between  users  both  at  30  deg  latitude  and 
separated  by  45  deg  in  longitude  and  can  be  found  from  the  plot  in 
Fig.  4.9a.  The  design  value  (kga)^^  can  then  be  read  from  the  plot 
in  Fig.  2.7. 

The  design  value  (X/Xq)^^  should  be  selected  to  typify  the 
average  value  that  might  be  experienced  by  user  pairs  located  anywhere 
within  the  service  area.  Failing  an  exhaustive  search  over  possible 
user  pair  locations,  an  adequate  approximation  can  be  obtained  from 
Figs.  4.10a  to  4.10c  by  noting  that  there  is  little  variation  during 
the  course  of  a  satellite  transit  for  a  given  user  pair  and  that  the 
results  are  not  sensitive  to  the  specific  user  pair  locations  con¬ 
sidered.  Therefore,  the  average  value  of  X/XQ  experienced  by  the  user 
pair  located  at  the  opposite  extremities  of  the  service  area,  as  given 
in  Fig.  4.10b,  will  be  used  for  the  design  value  (xA(pdes’ 

The  design  values  that  result  from  this  procedure  are  listed  in 
Table  4.1  for  reference.  The  orbital  altitude  of  3000  km  is  not  con¬ 
sidered  because  it  does  not  allow  coverage  of  the  entire  service  area. 

The  design  values  from  Table  4.1  and  the  operating  wavelength 
ratio  (X/XQ)  from  Eq.  (3.42)  can  be  used  to  determine  kga  from  Eq. 

(4.14).  Then,  using  Eq.  (4.9)  for  the  bistatic  scattering  angle  6  and 

2  2 
Eq.  (A.l)  for  o^/ita  »  the  normalized  scattering  cross  section  o/L  is 

calculated  from  Eq.  (4.13).  The  :sults  are  plotted  in  Figs.  4.11a  to 

4.1ic  for  the  user  pairs  at  the  extremities  of  the  service  area.  The 

altitude  dependence  stems  in  part  from  the  smaller  bistatic  scattering 
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Table  4.1 


TYPICAL  GENERALIZED  PACSAT  DESIGN  PARAMETERS 


Orbital 
Altitude 
h,  km 

Maximum  Bis tat ic 
Scattering  Angle 
6,  deg 

Design  Normalized 
Sphere  Size 

(k0a>des 

Design  Operating 
Wavelength  Ratio 

4000 

37.86 

1.1523 

0.7980 

5000 

33.43 

1.1412 

0.8417 

6000 

29.80 

1.1340 

0.8718 

7000 

26.81 

1.1287 

0.8936 

8000 

24.32 

1.1250 

0.9103 

9000 

22.23 

1.1221 

0.9233 

10,000 

20.45 

1.1200 

0.9337 

Satellite  longitude  relative  to  mid-user  longitude 


Fig.  4.1  la  —  Normalized  scattering  cross  section  for  a  nonlibrating  PACSAT 
in  a  circular  equatorial  orbit  as  viewed  by  users 
both  at  30  deg  latitude  and  separated  by 
45  deg  in  longitude 


Normalized  scattering  <5'  Normalized  scattering 

cross  section  cx/L2  cross  section  cr/L2 
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angles  but  mostly  from  the  larger  operating  wavelength  ratios  experi¬ 
enced  at  the  higher  orbital  altitudes.  The  variation  with  user  pair 
location  and  relative  satellite  longitude  is  small  at  a  given  alti¬ 
tude.  The  approximation 


o 


=  £=0. 6L2 


(4.15) 


is  probably  adequate  for  rough  estimates. 

To  illustrate  a  specific  design,  consider  a  system  designed  to 
operate  in  the  7900  to  8400  MHz  band.  Although  this  band  is  allocated 
for  earth-to-space  (i.e.,  uplink)  use  in  the  FIXED-SATELLITE  service 
(i.e.,  between  fixed  terrestrial  stations  and  an  active  satellite),  it 
would  be  suitable  for  use  by  a  passive  satellite  system  on  a  noninter¬ 
fering  basis.  The  wavelength  =  3.681  cm  corresponding  to  the 

midband  frequency  8150  MHz  could  then  be  used  to  determine  the  spe¬ 
cific  sphere  size  and  spacing  from  the  relationships 


a 


*mid 

2it 


(k0a)des' 


*mid 

2(171 


O^des 


(4.16) 


The  values  for  sphere  diameter  d  and  spacing  s  corresponding  to  the 
values  given  in  Table  4.1  are  listed  in  Table  4.2. 


Table  4.2 

TYPICAL  SPECIFIC  PACSAT  DESIGN  PARAMETERS 


Orbital  Altitude 
h,  km 

Sphere  Diameter 
d,  cm 

Sphere  Spacing 
s ,  cm 

4000 

1.351 

2.306 

5000 

1.337 

2.186 

6000 

1.329 

2.111 

7000 

1.322 

2.060 

8000 

1.318 

2.022 

9000 

1.315 

1.993 

10,000 

1.312 

1.971 
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The  slow  variation  of  sphere  diameter  with  orbital  altitude  re¬ 
flects  the  similarly  slow  variation  of  the  optimum  normalized  sphere 
size  with  bistatic  scattering  angle  displayed  in  Fig.  2.7.  A  depar¬ 
ture  from  the  listed  value  would  cause  the  reduction  in  normalized 
scattering  cross  section  associated  with  operating  away  from  the  area 
of  the  ridge  in  Fig.  2.6  and  might  therefore  not  be  significant.  A 
departure  from  the  listed  sphere  spacing,  however,  would  cause  a 
change  in  the  operating  wavelength  and  that  might  have  a  significant 
Impact  on  the  ability  of  the  system  to  remain  within  the  allotted 
frequency  band. 

Though  not  pertinent  to  the  topic  of  array  design,  it  is  conve¬ 
nient  to  list  in  Table  4.3  some  of  the  interesting  system  parameters 
associated  with  the  service  area  shown  in  Fig.  4.5.  The  minimum  arc 
of  mutual  visibility  is  found  from  Fig.  4.7c  and  the  minimum  mutual 
viewing  time  from  Fig.  4.8c  for  the  user  pair  at  the  upper  extremes  of 
the  service  area;  that  is,  both  users  at  45  deg  latitude  and  separated 
by  45  deg  in  longitude.  The  maxima  are  found  from  Figs.  4.7a  and  4.8a 
and  for  adjacent  users  at  the  lower  edge  of  the  service  area;  that  is, 
both  users  at  30  deg  latitude  with  no  separation  in  longitude. 


Table  4.3 

TYPICAL  PACSAT  SYSTEM  PARAMETERS 


Orbital 
Altitude 
h,  km 

Arc  of  Mutual 
Visibility,  deg 

Mutual  Viewing 
Time,  min 

Min.  No.  of 
Satellites 
Required 

Min. 

Max. 

Min. 

Max. 

4000 

14.3 

89.6 

7.9 

49.7 

26 

5000 

30.1 

99.3 

19.6 

64.6 

12 

6000 

41.4 

107.0 

31.3 

80.7 

9 

7000 

50.2 

113.2 

43.6 

98.3 

8 

8000 

57.3 

118.4 

56.8 

117.4 

7 

9000 

63.2 

122.8 

71.2 

138.4 

6 

10,000 

68.2 

126.6 

86.9 

161.3 

6 

The  minimum  number  of  satellites  required  is  the  integer  next 
larger  than  360  divided  by  the  minimum  arc  of  mutual  visibility.  This 
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is  equal  to  the  number  of  satellites,  which,  if  equally  spaced  about 
the  orbital  arc,  will  ensure  that  there  will  always  be  a  satellite  in 
mutual  view  of  all  users  within  the  service  area.  Inasmuch  as  PACSAT 
will  probably  not  have  a  station-keeping  ability,  the  satellites  will 
move  with  respect  to  one  another  because  of  initial  errors  in  orbital 
injection  velocity.  Thus,  a  practical  system  will  require  more  than 
the  minimum  number  listed  in  Table  4.3  if  continuity  of  service  is  to 
be  assured.  Anticipated  drift  rates  and  suitable  replenishment  strat¬ 
egies  for  preventing  incipient  gaps  will  be  considered  in  a  companion 
report . 

COMMUNICATION  CHARACTERISTICS 


Llne-of-Sight  Path  Length 

The  length  of  the  line  of  sight  from  the  user  to  the  center  of 
the  array  is  given  in  Eq.  (3.18).  The  maximum  value  of  the  central 
angle  0  is  given  by  Eq.  (3.17);  for  an  equatorial  orbit,  9  is  related 
to  the  satellite  longitude  and  the  user  latitude  and  longitude  by  Eq. 
(4.4),  as  illustrated  in  Fig.  4.1.  The  line-of-sight  path  length 
is  given  in  the  composite  plot  Fig.  4.12  as  a  function  of  these  param¬ 
eters. 

The  use  of  Fig.  4.12  is  illustrated  by  the  dashed  lines.  The 
plot  is  entered  at  the  upper  right  at  the  difference  beteen  the  satel¬ 
lite  and  user  longitudes.  At  the  intersection  with  tne  user  latitude, 
the  abscissa  is  followed  downward  to  the  intersection  with  the  satel¬ 
lite  altitude.  The  path  length  is  then  given  by  the  lower-left-hand 
ordinate.  As  noted,  the  abscissa  is  the  central  angle  6.  Here,  and 
in  the  following,  the  line-of-sight  path  length  is  designated  by  R 
rather  than  1^. 

The  Link  Equation 

In  Fig.  4.1,  let  the  users  be  designated  by  the  subscripts  T  and 
R,  rather  than  1  and  2,  to  denote  that  they  are  transmitting  and  re¬ 
ceiving,  respectively.  Let  P  denote  the  transmitter  power  and  G  the 
gain  of  the  antenna  each  is  using.  Then,  the  power  transmitted  in  the 


Line-of-sight  path  length  R,  km 
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direction  of  the  satellite  is  given  by  PTGT  and  the  power  density  at 

2  1 

the  satellite  by  ptgt/4itR^.  By  definition,  the  array  scatters  a 
power  pTGT°/4irR^  toward  the  receiver  with  unity  again.  Thus,  the 
power  density  at  the  receiver  is  given  by  PTG.jfl /(4it  )2r^r2  .  Let 
denote  the  effective  capture  cross  section  of  the  receiving  antennat 
The  received  power  is  then  given  by 


_ _ £« 

R  s2n2  _2 


(4*)^  Rj 


(4.17) 


where  L  is  a  loss  factor  (L  >  1)  that  accounts  for  pointing  inaccura¬ 
cies,  circuitry,  and  path  losses,  etc.,  in  practical  systems. 

Let  R  denote  the  signaling  rate  in  bits  per  sec.  Then,  the  bit 
duration  is  equal  to  1/R  and  the  received  energy  per  bit  Is  given 
by  Pr/R.  The  system  noise  power  spectral  density  NQ  is  given  by 

N0  -  kT  (4.18) 

where  k  ■  1.38054  x  10  2"*  Joules  per  degree  Kelvin  and  T  is  the  system 
noise  temperature.  The  received  signal-to-noise  ratio  is  then  given  by 


(4.19) 


Equation  (4.19)  can  be  written 


P-  >  P 


R 


min 


(4.20) 


which  states  that  the  received  power  PR  must  equal  or  exceed  a  mini¬ 
mum  detectable  value  PR  that  is  related  to  the  minimum  usable  signal- 

to-noise  ratio  (E./N.)  ml? 

min 
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The  factor  M  by  which  PR  exceeds  P  Is  called  the  system  margin. 

R  R  » 

Introducing  M  and  equating  Eqs.  (4.17)  anan(4.20)  then  leads  to 


M 


_ eff 

2  2  2 
(4ii)ZR^r‘l 


(4.21) 


which  characterizes  the  performance  of  a  passive  communication  link. 
The  effective  area  of  an  antenna  Is  related  to  Its  gain  by  the  rela¬ 
tionship 


G 


4x  A 


eff 


(4.22) 


where  \  Is  the  operating  wavelength.  Thus  Eq.  (4.21)  can  also  be 
written 


(4.23) 


For  a  circular  parabolic  dish  antenna  of  diameter  D  the  effective  area 
Is  given  by 


A 


eff 


n 


(4.24) 


where  n  Is  an  efficiency  factor  often  taken  as  0.55  for  practical 
dishes.  Then,  Eq.  (4.23)  can  be  specialized  to 


*  nT'18PTDX° 

64  X2R^R2L 


(A. 25) 
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>ical  Power  Budgets 


The  large  number  of  variables  present  in  Eqs.  (4.23)  or  (4.25) 
preclude  the  use  of  simple  plots  to  illustrate  system  performance  over 
a  wide  range  of  the  parameters.  The  power  budgets  listed  in  Table  4.4 
are  therefore  presented  to  show  what  results  from  using  typical  or 
practical  values. 


Table  4.4 

TYPICAL  PACSAT  LINK  PARAMETERS 
(PACSAT  in  circular  equatorial  orbit  at 
6000  km  altitude) 


Operating  frequency 

8  GHz 

8  GHz 

Transmitter  power 

5  kW 

10  kW 

Data  rate 

100  bps 

1  kbps 

Transmitter  antenna  diameter 

75  cm 

1  m 

Range  to  transmitter 

Scattering  cross  section 

10.000  km 
10®  m2 

10.000 
10®  m^ 

km 

Range  to  receiver 

10,000  km 

10,000 

km 

System  noise  temperature 

100  °K 

100  °K 

Losses,  margin,  and  threshold 

10  dB 

10  dB 

Receiver  antenna  diameter 

70  cm 

1.2  m 

System  operation  was  assumed  to  be  in  the  7.9  to  8.4  GHz  band 
cited  earlier.  Data  rates  of  100  bps  and  1  kbps  were  chosen  as  being 
representative  of  the  data  rates  of  interest  in  the  military  applica¬ 
tion  to  which  PACSAT  is  suited.  The  path  lengths  were  taken  arbitrar¬ 
ily  as  10,000  km;  the  path  length  to  the  horizon  is  10,600  km  for  an 
orbital  altitude  of  6000  km  (see  Fig.  4.12).  For  convenience,  the 

fi  7 

PACSAT  scattering  cross  section  was  taken  as  10°  m  ,  which  corresponds 

to  that  for  an  array  about  1300  m  long  (see  Eq.  (4.15)).  The  losses, 

margin,  and  threshold  LM  (E.  /N_)  .  were  lumped  into  a  single  factor, 

d  u  min 

taken  (perhaps  optimistically)  as  10  dB.  The  encouraging  aspect  of 
the  power  budgets  listed  in  Table  4.1  is  that  they  indicate  that  the 
link  can  be  supported  using  reasonable  transmitter  powers  and  antenna 
sizes.  Components  such  as  these  can  readily  be  accommodated  on  large 
aircraft  and  mobile  ground  terminals. 
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Appendix 

FIFTH-ORDER  R.M.S.  POLYNOMIAL  FIT  TO  MIE  SOLUTION 
EXPANSION  FOR  ELECTROMAGNETIC  SCATTERING  FROM 
A  PERFECTLY  CONDUCTING  SPHERE 


Mr.  V.  V.  Liepa  of  the  Radiation  Laboratory  at  the  University  of 

Michigan  at  Ann  Arbor  provided  calculations  of  the  real  and  imaginary 

2 

parts  of  Sj(0)  and  $2(6).  These  were  used  to  calculate  o^/x a  from 
Eq.  (2.18)  for  kga  *  0.9(0. 1)1. A  and  0  ■  0(10)50  deg.  A  uniformly 
weighted,  least-mean-square,  fifth-order  polynomial  in  kQa  and  6  was 
fitted  to  these  points  to  facilitate  computation  at  intermediate  points. 
The  polynomial  is  given  by 

8^/ira2  -  81.376  862  28  -  6.440  739  782  8  -  417.915  971  2  l^a 

+  14.772  132  26  62  +  26.659  069  49  9  kQa  +  814.742  881  3  (k0a>2 
+  0.361  421  426  2  63  -  42.417  652  88  9 2  kQa  -  40.368  661  59  6 (kQa)2 

-  741.250  998  2  (kQa)3  +  0.918  550  083  2  9  4  -  0.832  943  049  2  8  3kQa  (A.l) 
+  35.335  055  63  82(k0«)2  +  26.552  039  64  9(k0a)'*  +  320.020  340  2  (k0a)4 

+  0.040  049  287  89  95  -  0.944  447  017  0  94  kQa  +  0.463  573  159  1  83(k0a)2 

-  8.716  149  189  92(k0a)3  -  6.410  860  956  9 (kpa)4  -  53.335  225  21  (kQa)5 

where  6  is  in  radians.  The  r.m.s.  error  of  the  polynomial  at  the  36 

A  -3 

fitting  points  is  5.38  x  10  ;  the  largest  error  is  -  1.08A  *  10  ,  the 

smallest,  1.93A  *  10  The  polynomial  was  used  directly  to  plot 
9 

oT  /ia  in  Fig.  2.2;  in  accordance  with  Eqs.  (2.38)  and  (2.A1),  it  was 

2 

simply  multiplied  by  kQa  to  plot  o/L  (a/X)  in  Fig.  2.6. 
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